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Attached  is  a  user's  manual   on  the  HPR  Part  II  research  study  entitled 
"Design  of  Reinforced  Embankment s ."  This  report  is  the  implementation  part  of 
Task  i*    of  the  approved  work  plan.   The  authors  of  the  report  are  D.  N.  Humphrey, 
V.  0.  McCarron,  W.  F.  chen,  and  myself. 

The  user's  manual  describes  the  operation  of  the  finite  elenent  prograr 
PS-NFAP  developed  in  this  research.   The  program  analyzes  plane  strain  problems 
such  as  highway  embankments  on  soft  foundations  using  a  cap  type  soil  plasticity 
behavior  model.   The  use  of  four  other  programs  necessary  for  the  complete 
analysis  is  also  described.   Nos .  (1)  and  (2),  CPCALC  and  CAP,  calculate  the 
required  cap  model  parameters  from  the  results  of  typical  soils  tests.   No.  (3) 
AUTOGEN  generates  the  input  file  for  PS-NFAP,  and  No.  (A)  NFKINX  optimizes  the 
node  numbering  of  the  finite  element  mesh  used  in  the  problem.   PS-NFAP  has  a 
number  of  features  specially  adapted  for  the  design  of  reinforced  and  unrein- 
forced  embankments  on  soft  foundations.  A  detailed  example  illustrating  how 
the  series  of  programs  are  used  is  also  given  in  the  manual.   The  report  is 
essential  to  highway  engineers  contemplating  using  PS-NFAP  for  the  analysis  of 
embankments  on  soft  foundations. 

Copies  of  the  report  will  be  submitted  to  the  IDOH  and  FHVA  for  their  review, 
My  co-authors  and  I  look  forward  to  receiving  their  comments  on  the  manual. 
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This  user's  m.inual  describes  the  operation  of  PS-NFAP,  a  finite  element 
program  developed  for  the  analysis  of  plane  strain  problems  with  the  cap  soil 
behavioral  model.   In  addition,  two  programs  (CPCALC  and  CAP)  that  assist  in 
calculating  the  cap  model  parameters  from  typical  soil  tests  results,  a 
preprocessor  program  (AUTOCEN)  that  generates  the  PS-NFAP  input  file,  and  an 
auxiliary  program  (NFM1NX)  that  optimizes  the  FE  mesh  node  numbering  are 
described.   PS-NFAP  has  special  features  for  analysis  of  reinforced  and  unrein- 
forccd  embankments  constructed  on  soft  ground.  A  detailed  example  to  illustrate 
how  the  series  of  programs  are  used  to  analyze  a  reinforced  embankment  is  given. 

The  programs  a^e  written  in  FORTRAN  IV  and  FOKTRAK77.   They  are  implemented 
on  two  computer  systems:   (1)  IBM  3063  mainframe  using  VS-FORTRAN;  and  (2)  the 
series  of  IBM  personal  computers  with  a  math  co-processor  chip  using  Ryan-McFar lane 
FORTRAN,  version  2 .0  {Ryan-McFarland ,  1985).   The  IBM  personal  computers  include 
the  IBM-PC,  IBM-XT,  IBM-AT,  or  compatible  computers. 
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FE  ANALYSIS  OF  PLANE  STRAIN  PROBLEMS  WITH 
PS-NFAP  AND  THE  CAP  SOIL  MODEL 

USER'S  MANUAL 


INTRODUCTION 

This  user's  manual  describes  operation  of  PS-NFAP,  a 
finite  element  program  for  analysis  of  plane  strain  problems 
with  the  cap  soil  behavior  model.   In  addition,  two  programs 
(CPCALC  end  CAP)  that  assist  in  calculating  the  cap  model 
parameters  from  typical  soil  tests  results,  a  preprocessor 
program  (AUTOGEN)  that  generates  the  PS-NFAP  input  file,  and 
an  auxiliary  program  (NFMINX)  that  optimizes  the  FE  mesh 
node  numbering  are  described.   PS-NFAP  has  special  features 
for  analysis  of  reinforced  and  unrelnforced  embankments  con- 
structed on  soft  ground  (McCarron,  1985;  Humphrey  and  Holtz. 
1986).   A  detailed  example  that  Illustrates  how  the  series 
of  programs  are  used  to  analyze  a  reinforced  embankment  Is 
given.   Details  of  the  cap  model  are   given  in  McCarron 
(1985)  and  Humphrey  and  Holtz  (1986).   Procedures  to  deter- 
mine the  cap  parameters  from  soil  test  results  are  given  In 
Humphrey  and  Holtz  (1986).   Only  aspects  of  the  model  that 
deal  specifically  with  program  operation  are  discussed 
herein.   All  the  programs  use  the  convention  that  compres- 
sive stresses  and  strains  are   negative.   In  several  in- 
stances the  stresses  are   expressed  In  terms  of  the  first 
Invariant  I!  and  the  second  Invariant  of  the  stress 
devlator  tensor  J2   where 

i;  -  o;  4  c2  +  c3 

2  2  2 

J2  *  (or°2)  *    lo2~°3)      *  l03~0ll 

and  where:      o.  «  major  principal  stress 

o-  «  intermediate  principal  stress 

o-  «  minor  principal  stress 

The  programs  are   written  In  FORTRAN  IV  and  FORTRAN77. 
They  are    Implemented  on  two  computer  systems:  (1)  IBM  3083 
mainframe  using  VS-FORTRAN;  and  (2)  the  series  of  IBM  per- 
sonal computers  with  a  math  co-processor  chip  using  Ryan- 
McFarland  FORTRAN,  version  2.0  (Ryan-McFarland,  1985).   The 
IBM  personal  computers  Include  the  IBM-PC.  IBM-XT,  IBM-AT. 
or  compatible  computers  and  are   collectively  referred  to  as 
IBM-PC  computers  in  this  manual. 
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An  overview  of  the  steps  In  the  analysis  procedure  Is 
given  In  the  next  section.   A  description  of  each  program 
and  input  Instructions  are  given  in  the  Appendices.   CPCALC 
and  CAP  are  described  fn  Appendices  A  and  B,  respectively. 
AUTOGEN  is  described  In  Appendix  C.   PS-NFAP  and  NFH1NX  ere 
described  In  Appendix  D.   Appendix  E  shows  how  the  programs 
•re  used  to  analyze  reinforced  embankments  constructed  on 
•oft  ground;  also,  example  Input  and  output  files  for  each 
program  ere  given. 

OVERVIEW  OF  ANALYSIS  PROCEDURE 

The  first  step  In  the  analysis  procedure  Is  to  calcu- 
late the  cap  parameters  from  soil  test  results  using  the 
procedures  given  In  Humphrey  and  Holtz  (1986).   Three  of  the 
critical  parameters  ere  calculated  using  CPCALC  as  described 
In  Appendix  A.   In  addition,  there  are  several  cap  parame- 
ters which  control  execution  of  CAP  and  PS-NFAP.   These  ere 
described  in  Appendices  B  end  D. 

The  next  step  Is  to  compare  calculated  response  to 
behavior  exhibited  by  representative  laboretory  test  sam- 
ples.  For  example,  the  calculated  and  observed  stress 
strain  curves  from  a  triexiel  test  could  be  compared.   This 
Is  done  using  the  CAP  program  described  In  Appendix  B.   If 
necessary,  the  cap  parameters  can  be  adjusted  to  obtain  e 
better  fit  using  recommendations  given  in  Humphrey  end  Holtz 
(1986) . 

Next  the  problem  geometry  Is  defined  end  the  finite 
element  mesh  is  selected  using  the  example  in  Appendix  E  for 
guidance.   Then  the  input  file  for  PS-NFAP  Is  prepared  using 
the  instructions  In  Appendix  D.   The  preprocessor  program 
AUTOGEN  Is  used  to  generate  much  of  the  node,  element,  and 
load  data  as  described  in  Appendix  C.   If  the  ana  lysis  will 
be  carried  out  on  an  IBM-PC  the  euxlllery  node  renumbering 
program  NFM1NX  Is  run.   The  nodes  are  renumbered  by  PS-NFAP 
In  the  mainframe  version. 

The  final  step  Is  to  run  the  analysis  using  PS-NFAP. 
The  input  file  should  be  carefully  checked  for  errors.   The 
program  output  should  be  examined  to  be  sure  the  calculated 
results  are  reasonable.   The  steps  In  the  analysis  procedure 
ere  summarized  In  Fig.  1  end  illustrated  with  the  example  In 
Appendix  E. 
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Fig.  I.   Steps  In  analysis  procedure 
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CPCALC 
PROGRAM  TO  CALCULATE  CAP  PARAMETERS 
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CPCALC  Is  a  program  to  compute  three  critical  cap 
parameters,  namely:  the  hardening  parameters,  Dx   and 
W/(a-b).  and  the  aspect  ratio.  R.   The  program  prompts  the 
user  for  the  soil  properties  shown  In  Table  A.I. 

The  program  uses  the  solution  procedure  given  In  Humphrey 
•nd  Holtz  (1986).   The  values  of  *f/*0«  Rmax*  snd  (xf/xo)max 

•re  also  computed.   The  calculated  parameters  are   defined 
•nd  their  correspondence  with  the  variable  names  used  by 
CPCALC  is  shown  in  Table  A. 2.   Example  Input  and  output  Is 
shown  1 n  Append  1 x  E . 

REFERENCE 
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Tab  1 e  A . 1 
Input  soli  properties  for  program  CPCALC 

Soil     Variable  Description 

property     name 

o         am     slope  of  Drucker-Prager  surface 

1  /2 

r/o'         acn     normalized  J  '   Intercept  of 

vo  Drucker-Prager  surface 

j,/2/o'       sJ2fn   normalized  shear  strength  ratio 
2f    vo 

K         rkc     coefficient  of  lateral  earth  pressure 
°  at  rest 

(a-b)/b     abb    soil  compressibility  where  a  and  b  ere 

slopes  of  virgin  compression  and 
unloading- re  loading  curves  on 
e  -ln(o')  plot  from  hydrostatic 
consolidation  test 

BSCSSSSSSSSESKeSSSSCSSCeSCSESeSSSSESSCCESSSCSSKSKBESBSBBSS 

Note:  The  ratios  acn  and  sJ2fn  are  negative  since 
compressive  stresses  are  negative 
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Table  A.  2 
Cap  parameters  calculated  by  program  CPCALC 

Cap     Variable  Description 

parameter    name 

dlmens lonl ess  hardening  parameter 

normalized  hardening  parameter 

cap  aspect  ratio 

ratio  of  failure  to  Initial  1J 
intercept  of  cap 

maximum  possible  R  for  given  Input  soil 
propert les 

xfxomx   maximum  possible  x,/x   for  given  Input 
sol  1  propert les 


D*o 

dxo 

W/(a-b) 

wab 

R 

r 

*f/xo 

xf  xo 

R 
max 

rmax 

(W, 

rnax 

xf  xom 

APPENDIX  B 

CAP 
PROGRAM  TO  CALCULATE  SOIL  RESPONSE 
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The  CAP  computer  program  uses  the  cap  model  shown  In 
Fig.  6.1  to  calculate  soli  response  due  to  a  set  of  applied 
stresses  or  strains.   Details  of  the  cap  model  are  given  in 
HcCarron  (1985),  Chen  and  Batadl  (1985),  and  Humphrey  and 
Holtz  (I9B6).   The  CAP  program  was  developed  by  HcCarron  and 
Chen  (1986).   Its  main  use  Is  to  compare  calculated  behavior 
to  observed  test  sample  behavior.   For  example,  the  stress- 
Strain  curve  obtained  from  a  trlaxlal  compression  test  can 
be  compared  to  the  curve  calculated  by  the  program  for  a 
simulated  sample  with  the  same  Initial  conditions  subjected 
to  the  same  stress  path.   If  necessary,  the  cap  parameters 
can  be  adjusted  to  obtain  a  better  fit  using  the  guidelines 
given  tn  Humphrey  and  Holtz  (1986).   The  program  can  model 
ax  1  symmetric,  plane  strain,  and  plane  stress  conditions. 
Material  nonl Ineartty  or  both  material  and  geometric  nonltn- 
ear I ty  can  be  represented.   The  user  specifies  either  the 
stress  path  and   the  program  calculates  the  resulting  strains 
and  pore  pressure  or  the  strain  path  and  the  program  calcu- 
lates the  resulting  stresses  and  pore  pressure.   The  results 
are  output  In  tabular  form.   Stresses  and  pore  pressure 
versus  axial  strain  and  stress  path  data  are   stored  in  two 
files  for  convenient  Input  Into  commercial  plotting 
programs. 

The  program  Is  written  In  FORTRAN77  and  consists  of  a 
main  program,  15  subroutines,  and  a  block  data  unit.   It  is 
Implemented  on  an  IBM-PC  using  Ryan-McFarland  FORTRAN, 
Version  2.0  (Ryan-McFar land,  1985).   Program  execution  Is 
controlled  by  the  main  program.   The  purpose  of  the  primary 
subroutines  Is  summarized  In  Table  B. 1  and  the  purpose  of 
the  cap  model  subroutines  Is  summarized  in  Table  B.2. 
Although  the  cap  model  subroutines  have  the  same  names  as  In 
PS-NFAP  there  are   slight  differences  between  them  and  sub- 
routines from  CAP  and  PS-NFAP  are   not  Interchangeable.   The 
program  flowchart  Is  shown  In  Fig.  B.2.   Input  Instructions 
are   given  In  the  first  section  of  Appendix  B.l  followed  by 
comments  to  the  Input  Instructions.   All  Input  Is  free 
format.   Recall  that  compressive  stresses  and  strains  are 
negative.   To  execute  the  program  type: 

OCAP    <FILE.IN    >FILE.0UT 

where  FILE. IN  Is  the  name  of  the  Input  file  and  FILE. OUT  is 
the  name  of  the  output  file.   An  example  Input  file  Is  shown 
In  Appendix  E. 

The  calculated  strains  become  unbounded  at  failure. 
This  causes  error  messages  to  be  displayed  on  the  screen  or 
In  the  output  file.   When  operated  on  a  IBM-PC*  the  solution 
may  be  terminated  at  this  point  with  the  * control -break'  key 
or  ft  may  be  allowed  to  continue  until  the  Iteration  limit 
Is  exceeded. 
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The  output  file  lists  the  control  parameters,  cap  pa- 
rameters, and  calculated  stress  and  strain  components.   In 
addition,  two  files  are   generated  containing  stress-strain 
behavior  for  use  as  Input  Into  commercial  plotting  programs. 
The  first  file  Is  named  *plot.l'  and  contains  «-_• 

iozz'°yy)nZo^oU    Au/ovo'  °zz/oyy  *nd  4" ?/ovo 
In  columns  I  through  5,  where: 

t   ■  axial  or  vertical  strain 
zz 

o        «  axial  or  vertical  stress 

o   ■  radial  or  horizontal  stress 
yy 

au   •  excess  pore  water  pressure 

J_   «  second  Invariant  of  stress  devlator  tensor 

o'   «  Initial  axial  or  vertical  effective  stress 
vo 

The  second  file  Is  named  *plot.2'  and  contains  e   .  p ' ,  q, 

1/2 
I',  and  J_    In  columns  1  through  5,  where: 

p'   ■  (ozz+°yy>/2 

q   ■  (azz-°yy)/2 

I!   =  first  invariant  =  o'   *  o'   4  o' 
I  xx     yy     zz 

o'   *  lateral  or  horizontal  stress 
xx 

Note  that  '  Indicates  effective  stress.   Stresses  and 
strains  In  the  two  plot  files  are   output  using  the  soil 
mechanics  convention  that  compressive  stresses  and  strains 
are   positive.   An  example  output  1s  shown  In  Appendix  E.l. 
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Subrout Ine 

■SKCCECCEC1 

IN1T 

EL2D10 

RESET 

MATSET 

SOLVER 

FUNCT2 

WRITES 


Tab  1 e  B . 1 
Function  of  principal  subroutines 

Funct Ion 

Initialize  parameters  and  read  end  print 
material  data 

Access  to  cap  model 

Stores  current  material  state  on  disk 

Sets  coefficients  in  equilibrium  equations  and 
load  vector 

Forward/backward  substitution  of  equilibrium 
equat  ions 

Find  Interpolation  functions  and  their 
der i  vat  i  ves 

Print  stresses  and  strains  to  output;  write 
material  response  to  disk  files  'plot.l'  and 
*plot.2'  for  Input  to  plotting  programs 
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Subrout i ne 

«EBBB  =  =  X  =  =: 

EL2D10 
1EP210 
EP2D10 
EPL210 

PRAGER 

MD2D10 

HAXMIN 
DEVPRT 


:=cx=xxxr=r 


Table  B.2 
Function  of  cap  model  subroutines 

Funct Ion 

:bssbbcbebesbbssbxssebbsssbss=ssebe: 

Access  to  ceo  model 
Initialize  material  history 
Driver  for  material  model 


Compute  material  response  for  given  strain 
increment 

Form  constitutive  relation  for  Drucker-Prager 
material 

Form  constitutive  relation  for  work  hardening 
cap 

Determine  principal  stresses 

Determine  stress  invariants  and  deviatoric 
components 


TENCT 


Tens i le  loading 
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Figure   B.I.      Cap  model. 


-  B7 


START 
CALL  INI T 


LOOP  OVER  . 
STRAIN  PATH 


\ 

El 

I 


CALL  EL2D10 

I 

CALL  WRITES- 


LOOP  OVER 
STRESS  PATH 

CALL  EL2D10 

I 
CALL  RESET 

I 

CALL  MATSET 

I 
ITERATION  LOOP 

I 

CALL  RESET 

I 
CALL  MATSET 

I 

CALL  SOLVER 

I 
CALL  EL2D10 

I 

CONVERGENCE 

V 
YES 

k 
-CALL  WRITES 


•NO 


Figure  B.2.   Program  flowchart 
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INPUT  INSTRUCTIONS  FOR  CAP 
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INPUT  INSTRUCTIONS 

>>>  ALL  INPUT  IS  FREE-FORMAT  <<< 
Card  I :   HEADING 

Variable  Continent-  Entry 

TITLE        (1) 


Enter  the  title  for  use  In  labeling  the 
output 


•  Consents  to  the  Input  Instructions  are   given  In  the 
f o 1  I ow I ng  sect  I on 


Card  2: 

Variable 

1TYP2D 


ANALYSIS   INFORMATION 


INDNL 


Comment 


(2) 


Entry 


(3) 


Analysis  type; 
.EQ.  0;  axlsymmetrlc  case 
.EQ.  I;  plane  strain  case 
.EQ.  2;  plane  stress  case 

Type  of  nonlinear  analysis; 
.EQ.  1;  material  nonllnearlty  only 
•EQ.  4;  both  material  and  geometric 
nonl I  near  I ty 
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Card  3: 

CAP  PARAMETERS 

Variable 

Comment 

PROP( 1) 

(4) 

Kl 

PR0P<2) 

(4) 

K2 

PR0P(3) 

(4) 

A 
P 

PR0P(4) 

(4) 

ml  i 

PROP ( 5 ) 

(5) 

Gl 

PROP ( 6 ) 

(5) 

G2 

PROP ( 7 ) 

(6) 

AM 

PROP ( B ) 

(6) 

AC 

PROP(9) 

(7) 

AW 

PROP (10) 

(7) 

AD 

PROP (11) 

(8) 

R 

PROP( 12) 

(9) 

XL 

PROP (13) 
PROP (14) 


PROP( 16) 
PROP( 17) 

PROP (18) 
PROPU9) 


(10) 
(11) 


PROP (15)     (12) 


(13) 
(14) 

(16) 
(15) 


Entry 

-  bulk  modulus  parameter 

-  bulk  modulus  parameter 

-  atmospheric  pressure 

-  minimum  bulk  modulus 

-  shear  modulus  parameter 

-  shear  modulus  parameter 

-  slope  of  Druckei — Prager  surface 

1/2 

-  J_    Intercept  of  Drucker- 

Prager  surface 

-  hardening  parameter 

-  hardening  parameter 

-  cap  aspect  ratio 

-  I'  coordinate  of  Intersection 
of  cap  and  Drucker-Prager 
surface 

.EQ.  "fl.i  normally  consolidated 
sol  1 

TENCUT  -  limiting  tensile  strength 

STATE   -  movement  of  cap 

.EQ.  It  expansion  and 

contraction 
.GT.  1;  expansion  only 

Al      -  unit  weight  of  soil; 
must  be  0. 

A2     -  Initial  vertical  stress 

K       -  ratio  of  Initial  horizontal  to 
vertical  stress 

open   -  not  used;  must  be  0. 

FAC    -  pore  pressure  response  factor; 
.EQ.  O.i  drained 
.GT.  O.t  undralned 
.LT.  0.;  undralned 
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Card  3:  CAP  PARAMETERS  (CONTINUED) 

Variable  Comment              Entry 

PROP(20)  (16)    T1MEON  -  not  used;  must  be  0. 

PROP(21)  (16)    WGT    -  not  used;  must  be  0. 

PROP (22)  (16)     BOUY    -  not  used;  must  be  0. 


Card  4:   SOLUTION  TYPE 
Va  r t  ab 1 e   Comment 


Entry 


SOLTYP 


(17)     Solution  type 

.EQ.  'strain' »  strain  path  specified 
.EQ.  %stress';  stress  path  specified 


FOR  STRAIN  PATH  SPECIFIED 
Card  5:   STRAIN  PATH  DATA 
Variable    Comment 
NTIMES       (18) 


DSTRAIN(l)  (IB) 

D5TRAIN(2)  (18) 

DSTRAINO)  (18) 

D5TRAIN(4)  (18) 

DSTRA1N(5)  (19) 


Entry 


Number  of  times  to  apply  strain 
Increment 

Ac   5   radial  or  horizontal  strain 
^    Increment 

Ae   ;   axial  or  vertical  strain 


zz 


Increment 


A2e   J  engineering  shear  strain 
**      Increment 

At   ;   radial  or  horizontal  strain 


xx 


Aw; 


Increment 

spin  Increment  for  large 
displacement  analysis 


Repeat  card  5  for  each  set  of  strain  Increments.   Input  Is 
terminated  by  a  card  5  with  NTIMES-0  and  DSTRAINd )-0. I  this 
must  be  the  last  card  In  the  Input  file. 
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FOR  STRESS  PATH  SPECIFIED 
Card  5:   ITERATION  INFORMATION 

Entry 


Variable 

Comment 

RTOL 

(20) 

1TEMAX 

(20) 

IREF 

(21) 

Solution  tolerance  for  convergence; 

■  |  !  unbalanced  stress,' J/J  {total  stress!! 

max  1 mum  number  of  Iterations  allowed 

maximum  number  of  reformations  of  the 
constitutive  relation 


Card  6:   STRESS  PATH  DATA 

Variable    Comment  Entry 

KSWTCH       (22)     Constitutive  relation  flag 

.EQ.  0;  cap  model  constitutive  relation 
.EQ.  1;  elastic  constitutive  relation 

NTIMES       (22)     Number  of  times  to  apply  stress 

Increment 

DSIG(l)      (22)     La       i    radial  or  horizontal  stress 

'y   Increment 

DS1G(2)      (22)     La      x    axial  or  vertical  stress  Increment 

DSIGO)      (22)     La      ;  shear  stress  Increment 

yz 

DS1G(4)      (22)     La       ;  radial  or  horizontal  stress 

xx   Increment 

Repeat  card  6  for  each  set  of  stress  Increments.   Input  Is 
terminated  by  a  card  6  with  NTIMES=0  and  DS1G(I)-0.1  this 
must  be  the  last  card  In  the  Input  file. 
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COMMENTS  FOR  CAP  INPUT  INSTRUCTIONS 

1.  The  TITLE  Is  used  to  label  the  program  output.   It  may 
be  up  to  80  characters  long. 

2.  ITYP2D  .EQ.  0  Is  used  for  trlaxlel  test  conditions.   For 
this  case  the  Z-axIs  Is  the  axial  direction  and  the 
Y-axis  Is  the  radial  direction.   ITYP2D  .EQ.  1  Is  used 
for  plane  strain  test  conditions.   For  this  case  no 
strain  Is  allowed  In  X-di recti on;  the  Z-axts  Is  the 
vertical  direction  and  the  Y-axis  Is  the  horizontal 
direction.   For  1TYPE2D  .EQ.  I,  stress  In  the 
X-dlrectlon  Is  constant. 

3.  INDNL  should  be  the  same  as  will  be  used  In  PS-NFAP. 
For  most  cases  both  material  and  geometric  nonllnearlty 
will  be  considered  so  INDNL  .EQ.  4. 

4.  The  bulk  modulus  (K)  is  computed  according  to  the 
f o 1 1 ow 1 ng  equat 1 on : 

K  .  HAX(K|*Ap-I-i;/fS-ApnK2-,  *nin> 

K  and  K_  are  fitting  parameters.  A   Is  atmospheric 

pressure  and  Is  used  to  obtain  a  dlmens lonless  equation, 

and  K  ,   Is  the  minimum  value  of  K  which  will  control 
mln 

for  small  I'.      K  ,   must  be  .GE.  0. 
1     mi  n 

5.  The  shear  modulus  (G)  Is  computed  according  to  the 
f o 1 1 ow I ng  equat I on : 

G  «  Gj  ♦  G2«K 

G.  and  G?  are   fitting  parameters. 

6.  AM  and  AC  are  defined  in  Fig.  B.l. 

7.  AW  and  AD  are   material  parameters  for  the  hardening  rule 
which  is  given  by  the  following  equation: 

ep  «  AW*(exp(AD-X)  -  I.] 

where  X  Is  the  Intersection  of  the  cap  with  the  Ij-axls. 
For  large  X,  ep  approaches  AW. 

8.  R  Is  the  aspect  ratio  of  the  work-hardening  elliptical 
cap.   R  Is  defined  as  the  ratio  of  the  horizontal  to 
vertical  axis  of  the  ellipse. 

9.  XL  .EQ.  +1.  for  normally  consolidated  soil  and  the 
program  calculates  the  Initial  cap  position. 
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10.  TENCUT  Is  the  maximum  principal  tensile  strength.   The 
material  is  assumed  to  fracture  after  reaching  this 
value.   For  most  test  conditions  In  soli  mechanics,  the 
principal  stresses  remain  compressive  so  the  value  of 
TENCUT  Is  unimportant  and  It  may  be  set  .EQ.  ♦  !. 

11.  STATE  Is  a  material  parameter  which  controls  whether  the 
cap  Is  allowed  to  contract  as  well  as  expand. 

12.  Al  Is  the  unit  weight  of  soil.   It  must  be  0  for  use  In 

CAP. 

13.  A2  Is  used  to  specify  the  initial  vertical  (Z-dlrect Ion) 

stress  (o'  ). 
vo 

14.  K   Is  used  to  compute  the  Initial  horizontal  (Y  and  X- 
d°rectlons)  stress  (<»I)  according  to  the  following 

f ormu I a : 

o'   «  K   •  a' 

ho     o     vo 

For  hydrostatic  (Isotropic)  Initial  conditions  K   ■  1.0. 

o 

15.  FAC  Is  used  to  compute  the  bulk  modulus  of  water  (K  ) 
for  use  In  undralned  analysis:  w 

K   «  FAC  •  K 

w 

FAC  .EQ.  0.  for  drained  analysis 

.GT.  0.  for  undralned  analysis;  computes  both 

positive  and  negative  excess  pore  pressures 
.LT.  0.  for  undralned  analysis;  computes  only 

negative  (compressive)  excess  pore 

pressures;  positive  pore  pressures  set 

.EQ.  0. 

FAC  .EQ.  10.  has  been  found  suitable  for  undralned 
analysts  of  soil  mechanics  problems. 

16.  PROP(IB),  PROP(20),  PR0P(21),  and  PROP(22)  are   not  used 
by  CAP  and  must  be  0. 

17.  The  strain  path  may  be  specified  and  CAP  will  compute 
the  resulting  stresses  or  the  stress  path  may  be  speci- 
fied and  CAP  will  compute  the  resulting  strains.   The 
format  of  the  remaining  cards  depends  on  the  option 
chosen.   *straln'  or  *stress'  should  be  Input  In  lower 
case  without  the  quotes. 


-  B14  - 

18.  This  card  Is  used  If  the  strain  path  Is  specified.   It 
should  be  repeated  for  each  set  of  applied  strain  Incre- 
ments.  Smaller  increments  should  be  used  as  failure  is 
approached.   It  may  be  necessary  to  adjust  the  strain 
Increments  based  on  the  approximate  failure  strain  Indi- 
cated by  a  preliminary  run. 

For  axlsymmetric  conditions  (1TYP2D  .EQ.  0)  At    » 

yy 

ac   .   For  plane  strain  conditions  (1TYP2D  .EQ.  1) 

e   «  Ac    «  0. 
xx      xx 

19.  w  accounts  for  rigid  body  motion  (McCerron, 1985) ;  1t  Is 
not  accounted  for  In  most  soil  mechanics  tests  and  may 
be  taken  as  0. 

20.  This  card  and  the  following  card  are  used  If  the  stress 
path  Is  specified.   RTOL  .EQ.  0.005  and  ITEMAX  .EQ.  20 
are  suitable  for  many  test  conditions  encountered  In 
sot  I  mechanics. 

21.  IREF  specifies  the  number  of  times  the  stiffness  matrix 
will  be  reformed  If  the  solution  does  not  converge  In 
ITEMAX  Iterations.   If  the  solution  still  does  not 
converge  after  IREF  reformations  execution  Is 
terminated. 

22.  This  card  should  be  repeated  for  each  set  of  applied 
stress  Increments.   Smaller  Increments  should  be  used  as 
failure  Is  approached.   It  may  be  necessary  to  adjust 
the  stress  Increments  based  on  the  approximate  failure 
stress  Indicated  by  a  preliminary  run.   In  general 
KSWTCH  .EQ.  0  except  for  first  Increment  of  an  unloading 
sequence  that  would  cause  movement  of  the  state  of 
stress  from  the  cap  or  Druckei — Prager  surface  Into  the 
elastic  region  In  which  case  KSWTCH  .EQ.  1. 

For  axlsymmetric  conditions  (ITYP2D  .EQ.  0)  Ao    ■ 

Ao   .   For  plane  strain  conditions  (ITYP2D  .EQ.  1) 
xx 

La         ■  0. 
xx 
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AUTOGEN 
PREPROCESSOR    FOR    PS-NFAP 
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AUTOGEN 
PREPROCESSOR  FOR  PS-NFAP 

by  W.  0.  HcCarron  and  W.  F.  Chen 
with  revisions  by  D.  N.  Hunphrey 


AUTOGEN  commands  are  separated  into  three  categories:  JOINT, 
ELEMENT,  and  LOAD  commands.  The  global  categories  should  appear 
In  the  order: 

JOINT 

ELEMENT 

LOAD 


Two  auxiliary  commands  are  available  vhlch  allow  cards  sot 
generated  by  AUTOGEN  to  be  Incorporated  at  the  appropriate  place 
In  NFINP.   The  command  NFAP  causes  the  Immediately  following 
cards  to  be  copied  directly  into  NFINP  as  they  appear  and  In  the 
position  they  appear.   The  cards  are  copied  to  NFINP  until  a 
JOINT,  ELEMENT,  or  LOAD  command  is  encountered.   The  command  FOLL 
has  the  same  function  as  NFAP,  except  that  It  la  used  to  append 
the  last  group  of  cards  to  the  end  of  NFINP.   A  FOLL  command  and 
any  cards  to  be  appended  to  NFINP  must  be  the  final  card  sequence 
Id  the  AUTOGEN  INPUT  file. 

Input  Instructions  are  given  In  Appendix  C.l.   All  Input  Is 
formatted.   Sample  Inut  and  output  Is  given  in  Appendix  E. 


OPERATION 

READ  FIRST  CARD 

EXECUTE  JOINT  COMMANDS 

EXECUTE  ELEMENT  COMMANDS 

EXECUTE  LOAD  COMMANDS 

COMPLETE  KFAP  DECK 

STOP 
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SUBROUTINE 
MAIN 

JOINT  

ELEMNT  

LDCEN   DER1Q 

I 
OUTPT        FUNCT2 

I 
MAIN  OUTPT 


SUBkOUNTlNES    CALLED 


•ELCEN 

I 
ELDUP 

I 
ELSTRC 

I 
EMB2 

I 
REDEF 

I 
OUTPT 


EMB2 

I 

OUT1 


Fig.     C.l       AUTOCEN    Flowchart. 
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INPUT  INSTRUCTIONS  FOR  AUTOGEN 
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FIRST  CARD 

INPUT    :  numjt,  numel 
FORMAT   :  (215) 

Definitions : 

nun J t  -  Number  of  joints  In  the  model. 

Dunel  •  Number  of  two-dimensional  elements  In  the  model. 

The  remainder  of  the  cards  may  be  JOINT,  ELEMENT,  LOAD,  NFAP,  or  FOLL 
cards  as  described  below. 


EXAMPLE 

15     8 
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JOINT  COMMANDS 

The  Joint  definition  cards  and  joint  generation  comands 
should  be  preceded  by  a  JOINT  header  card.   Any  number  of  JOINT 
header  cards  aay  be  present  within  the  INPUT  file.   The  various 
options  available  for  joint  generation  are: 

JTCEN   :  To  generate  a  ctrlng  of  nodes  from  a  previously 
defined  anchor  node. 

JTDUP   :  To  create  a  new  set  of  nodes  from  a  previously 
defined  set  of  nodes. 

MSHGEN  :  To  generate  a  mesh  of  nodes  beginning  with  a 
previously  defined  node. 


EMB 


To  generate  a  nesh  of  nodes  for  an  embankment 
nodel.   Symmetry  about  the  centerllne  is  assumed, 
therefore,  only  half  of  the  embankment  is  generated. 
A  2  to  1  (h/v)  side  slope  of  the  embankment  Is 
assumed  . 


JTP1X   :  To  specify  the  type  of  restraint  for  a  group  of 
nodes. 


In  some  cases  it  may  be  convenient  to  redefine  a  node  defined  by 
a  previous  command.   The  most  recent  definition  Is  retained. 
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JOINT  HEADER   and  JOINT  DEFINITION  CARDS 

JOINT  DEFINITION  CARDS  :  To  define  a  nodal  point 

The  JOINT  header  card  say  be  followed  by  any  number  of 
JOINT  DEFINITION  cards. 

INPUT   :  node,  x,  y,  z 
FORMAT  :  (lx, K,,3f 10.0) 


Definitions 

node 

x 

y 

z 


•  The  node  number. 

*  X  coordinate  (leave  blank). 
■  Y  coordinate. 

-  Z  coordinate. 


EXAMPLE 

Joint 
1 
6 


-2. 
-6. 


0. 
3. 
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JTGEN 


To    generate    a    string    of    nodes    from    •    previously 
defined    anchor    code. 


INPUT 
FORMAT 


Def lni  t Ions 

en 

numy 

nlncr 


:  nn ,  nuny  ,  nlncr,  d 
:  (lx,14,215,2f 10.0) 


dy  ,  dz 


<Jy 

dz 


A  previously  defined  node  at  vhlch  the 

node  string  will  originate. 

The  total  number  of  nodes  In  the  string, 

Including  the  original  node. 

The  nodal  lncrenent  used  to  obtain 

node  numbe  rs . 

The  Y  coordinate  Increment  between 

The  Z  coordinate  Increment  between 


the  new 

nodes  . 
nodes  . 


EXAMPLE 


Jtgen 

1 


-1 
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JTDUP 


To  create  a  new  aet  of  nodes  from  a  previously 
defined  aet  of  nodes. 


INPUT    :  nn,  o2,  nlncr,  dy , 
FORMAT   :  ( lx , 14 , 21 5 , 2f 1 0.0) 


dz 


Def inlt ions 


nn 
d2 
nlncr 

dy 
dz 


•  The  first  node  of  the  previously  defined  aet 
■  The  last  node  of  the  previously  defined  aet. 

•  The  nodal  Increment  used  to  obtain  the  new 


node  numbers. 

The  Y  coordinate  increment. 

The  Z  coordinate  increment. 


EXAMPLE 

J  tdup 
1 
4 
7 

10 


3 

6 

9 

12 


3 
3 
3 
3 
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MSHGN   :  To  generate  a  mesh  of  nodes  beginning  with  • 
previously  defined  node. 

INPUT    :  nn,  numy ,  nyi,  dy ,  numz ,  nzl,  dz 
FORMAT   :  (lx, 14, 2(215, f 10.0)) 

Definitions : 

nn  -  The  previously  defined  anchor  node. 

numy  •  The  total  number  of  nodes  In  the  Y  direction, 

nyl  -  The  node  number  increment  In  the  Y  direction, 

dy  «*  The  coordinate   Increment  in  the  Y  direction. 

nunz  -  The  total  number  of  nodes  in  the  Z  direction, 

nzl  -  The  node  number  Increment  In  the  Z  direction, 

dz  -  The  coordinate   Increment  In  the  Z  direction. 


EXAMPLE 


■  shg 

1 

10 

U3 

1&B 


9 
8 

3 


14 

U 
8 
6 


6 

6 

12 

12 


9 

1 

5 

1 

5 

1 

3 

1 

-3 
-6 

-6 
-12 
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EMB 


To  generate  a  aeth  of  nodes  for  an  embankment 
model.   Symmetry  about  the  centerllne  Is  assumed, 
therefore,  only  half  of  the  embankment  is  generated 
A  2  to  1  (h/v)  aide  slope  of  the  embankment  Is 
a  ssumed. 


INPUT    :  nn,nel,nlyrs,dy,elv 
FORMAT   :  (315,  2fl0.0) 


Definitions : 

on    •  The  number  of  the  first  node  to  be  generated, 
nel   -  The  number  of  elements  in  the  first  row  (default 
nlyrs*  The  number  of  element  layers  in  the  embankment, 
dy    ■  Half  the  length  of  one  element, 
elv   -  The  elevation  of  the  first  node. 
•  Elevation  of  top  of  embankment. 
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EXAMPLE 
emb 

1     5 


NOTES 


(1) 
(2) 
(3) 

(4) 

(5) 

(6) 

(7) 


This  option  may  only  be  used  once  in  the  model. 

The  mesh  generated  is  for  eight-noded  elements. 

At  least  one  Joint  must  be  defined  prior 

to  this  command. 

Joint  restraints  in  the  y-dlrection  along  the  axis 

of  symmetry  are  automatically  generated. 

The  elements  have  dimensions  2dy  in  the  y-directlon 

dy  in  the  r-direction. 

first  node  (nn)  i6  located  on  the  centerllne  at 

top  of  the  embankment. 

usefulness  of  this  option  is  limited  if 
incremental  construction  of  the  embankment  will 
be  used  in  PS-NFAP  since  the  resulting  lift  thickness 
is  too  large  unless  an  excessively  large  number  of 
elements  are  used. 


and 
The 
the 
The 
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JTFIX   :  To  specify  the  type  of  restraint  for  *  group  of 
nodes 


INPUT   :  iope,  nfxy,  nfxz,  nd(i) 
FORMAT  :  (lx, 14,1515) 


Def init Ions : 


lope 


nfxy 
nfxz 

nd(i) 


0  -  Fix  all  nodes  between  consecutive  pairs. 

1  •  Fix  listed  nodes. 

2  "  Fix  nd(2)  number  of  nodes  beginning  with 

nd(l)  and  having  a  nodal  Increment  of  nd ( 3 ) 

•  Type  of  restraint  In  the  Y  direction. 

•  Type  of  restraint  in  the  Z  direction. 

-  The  list  of  nodes,  up  to  twelve  on  each  card 


EXAMPLES 


Jtf  ix 

0 
1 
2 
2 


0 
1 

0 

1 


1 

3 
16 
23 


3 

6 

19 

18 


13 
9 
3 
3 


15 
12 


15 
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ELEMENT  COMMANDS 

The  elenent  commands  are  used  to  generate  tvo-dlmens 1  one  1  4 
to  8  node  quadrilateral  elements.   The  element  definition  and 
generation  commands  should  be  preceded  by  an  ELEMENT  header  card. 
Any  number  of  ELEMENT  header  cards  may  be  present  within  the 
INPUT  file.   The  virlouE  options  available  for  element  generation 
•  re : 

ELDUP   :  To  create  a  new  group  of  elements  from  a  previously 
defined  group. 

ELGEN   :  To  create  an  element  mesh  from  an  existing  element. 

ELSTRG  :  To  create  a  string  of  elements  from  an  existing 
e  lenient  . 

EMB     :  To  generate  a  nesh  of  elements  for  an  embankment 
model.   Symmetry  about  the  centerllne  is  assumed, 
therefore,  only  half  of  the  embankment  is  generated. 
A  2  to  1  (h/v)  side  slope  of  the  embankment  is 
a  ssumed . 

REDEF   :  To  redefine  the  nodal  Incidences  of  a  string  of 
elements 
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ELEMENT  HEADER  ind  ELEMENT  DEFINITION  CARDS 

The  ELEMENT  header  card  may  be  followed  by  any  suaber 
of  ELEMENT  DEFINITION  cards. 

ELEMENT  DEFINITION  CARDS 

CARD  1    INPUT    :  nn  ,  lps  ,t)t  yp  ,bet  a  ,  thic 
FORMAT   :  (1  x.  14  ,215  ,2f  10.0) 


CARD  2 


INPUT 
FORMAT 


(nd(l).i-l ,8) 
(815) 


De  f  1  nl  t  1  ons 
CARD  1 


nn 
lps 
■  typ 
beta 

thlc 


-  The  element  number 

■  The  stress  output  table  (not  used;  enter  1) 

-  The  material  property  number  ^J 

-  The  orientation  angle  for  orthotroplc  tatrrlil 
(not  used;  leave  blank) 

■  The  element  thickness  for  plane  stress  eleme  r.^sf 
(not  used;  leave  blank) 


CARD  2    nd(l)    -  the  node  numbers  that  form  the  element 

If  node  position  not  used 

EXAMPLE 


eleme 
1 
1 


1 
2 


3 
5 


-  C13  - 


ELDUP   :  To  create  a  new  group 
defined  group. 


of  elements  from  a  previously 


INPUT 
FORMAT 


nn  , 
(lx, 


n2, 
14. 


nell  , 
315) 


njl 


Def  lnl t Ions : 

on    -  The  first  element  of  the  original  group. 
d2    -  The  last  element  of  the  original  group. 
nell  •  The  element  number  increment  used  to 

obtain  the  new  elements, 
njl   ■  The  node  number  Increment  used  to  obtain 

the  new  element  s . 


EXAMPLE 


eldup 

1 
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ELGEN 


To  create  an  element  mesh  from  an  existing  element. 


INPUT 
FORMAT 


nn,  nurny ,  nelly,  ndly,  nun:,  sells,  sdlz 
(  lx,  14,  615) 


Definitions : 


nn     -  The  previously  defined  anchor  element. 

sumy   -  The  total  number  of  elements  in  the  Y  direction 

including  anchor  element, 
nelly  •  The  Increment  in  element  number  in  the  Y  direction 
ndly   -  The  node  number  Increment  In  the  Y  direction, 
numz   -  The  total  number  of  elements  In  the  Z  direction 

Including  anchor  element, 
nellz  -  The  Increment  In  element  number  In  the  Z  direction 
ndlz   -  The  node  number  Increment  In  the  Z  direction. 


EXAMPLE 

e  lgen 
1 


-  CI5  - 


ELSTRC  :  To  create  •  string  of  eleoentE  from  an  existing 
e leoent  . 

INPUT    :  nn,  dud,  nell,  nd(i) 
FORMAT   :  (lx,  14,  1015) 

Def lnl t Ions : 

on    -    The  original  element. 

dub   ■  The  total  number  of  elenents  In  the  string  Including 

the  original . 
nell  •  The  Increment  In  the  element  number. 
nd(i)  -  An  array  of  node  number  Increments  to  be  added  to 

the  corresponding  node  number  of  the  element. 

EXAMPLE 


elstr 

1 
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EHB 


CARD  1 
CARD  2 


To  generate  a  aesh  of  elements  for  an  embankment 
sodel.   Symmetry  about  Che  centerllne  la  assumed, 
therefore,  only  half  of  the  embankment  la  generated 
A  2  to  1  (h/v)  aide  alope  of  the  embankaent  la 
aaauned . 

INPUT    :  nn , ips  .beta  .thick 
FORMAT   :  215,  2M0.0 

INPUT    :  (mt(l) ,1»1 ,nlyrs) 
FORMAT   :  1615 


Definitions : 


nn 


-  The  number  to  be  assigned  to  the  first  element 
lps   -  Stress  printing  option  (not  used;  enter  1). 
beta  -  Element  orientation  for  orthotroplc  material 

(not  used;  leave  blank), 
thick*  Element  thickness  If  plane  stress  option 

(not  used;  leave  blank), 
mt    -  The  material  type  to  be  assigned  to  each  layer 
nlyrs-  The  number  of  element  layers  In  the  embankment 


EXAMPLE 


emb 


NOTES 


1 
1 

(1) 
(2) 
(3) 


This  option  may  only  be  used  once  In  the  INPUT  file 

The  elements  generated  are  eight-noded  elements. 

A  least  one  element  must  be  defined  prior 

to  using  this  command  (a  dummy  element  may 

be  used  which  will  be  redefined  by  this  command). 
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REDEF   :  To  redefine  the  nodal  Incidences  of  a  string 
of  elements 

INPUT    :  nn,n2,nl ,(nd(l ) ,1-] ,8) 
FORMAT   :  Ix, 16, 215, 815 

Definitions: 


nn 

b2 
nl 

nd 


The  first  element  of  the  string. 

The  lsst   element  of  the  string. 

The  Increment  to  be  added  to  the  substituted 

nodes  nd(l ) . 

The  arrsy  of  nodes  to  be  substituted  for  element  "nn' 

(s  sero  entry  csuses  the  node's  original  definition 

to  remain  the  ssme). 


EXAMPLE 


redef 

147   158 


20 


21 


14 
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LOAD  COMMANDS 

AL'TOCEN  may  be  uaed  to  generate  gravity  loada  or  oqulval 
nodal  loada  to  balance  apeclfled  Initial  itrenci.   Naylor, 

•  a  a  ak*  a  i  liaaj    -    ~  _         •»«_.__.._  M  _  aat.  *.   _    _  * 


I 


ent 

t 

n- 
f  1- 

ng 
hor- 


CARD    1 
CARD    2 


INPUT 
FORMAT 
INPUT 
FORMAT 


LOAD 

(•4) 

ldtyp,ngrp,ncury,ncurz 

(4110) 


Definitions 


ldtyp  -  1  To  generate 

-  2  For  ln-ai  t  u 

-  3  For  ln-altu 

•  l*    For  In- situ 

•  5  To  generate 

nodes  In  an 


body  force  or  gravity  loads. 
loads  (y  direction  only), 
loads  (y  and  z  direction), 
loads  (z  direction  only), 
a  z-dlrectlon  unit  load  for  the 
element . 


ngrp 

o  cur  y 
n  c  ur  z 


The  number  of  element  groups 
(ngrp  nutr.ber  of  Card  3  must  be  provided). 
Load  curve  number  for  y-dlr  loads  (default 
Load  curve  number  for  z-dlr  loads  (default 


1) 
1) 


If  ldtyp  -  1 

CARD  3    INPUT    :   n ]  , n 2 , b f  ,  nnl nt z  ,  1 d 1 r 
FORMAT   :   ( 2 1 1 0 , f J  0 . 0 , 1 1 0  ) 

Definitions: 


ol     -  Flrat  element  number  of  this  group;  default 
n2     -  Last   element  number  of  this  group;  default 
bf     -  Body  force  (weight  of  aoil). 
nnlnt  -  Order  of  gauas  quadrature,  default  -  2. 
ldlr   -  Direction  of  load,  default  -  3. 

-  2   y  direction. 

-  3   c  direction. 


1  . 
nume  1 


If  ldtyp  -  2 

Same  n  ldtyp 
are  computed. 


■  3  except  only  z-dlrectlon  equivalent  loads 
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LOAD  COMMANDS  (CONTINUED) 
If  ldtyp  -  3 

CARD  3    INPUT    :   nl ,n2,al ,«2 ,xk  ,nnint 

FORMAT   :   (21 10, 3f 10.0,110) 

Definitions : 

nl     -  First  element  number  of  this  group;  default  ■  1. 
n2     -  La6t   element  number  of  this  group;  default  -  nuirel 
el,a2  •  Coefficients  used  to  compute  the  Initial  stresses, 
xk 

sz  -  al*z  ♦  a2 

ay  ■  s**xk 
nnlnt  -  Order  of  gauss  quadrature,  default  ■  2. 

If  ldtyp  -  4 

Same  86  ldtyp  •  3  except  only  r-directlon  equivalent  loads 
are  computed . 

If  ldtyp  -  5 

CARD  3    INPUT    :   nl ,n2 
FORMAT   :   (2110) 

Definitions : 

nl      •  First  element  number. 
n2      -  Last   element  number. 


EXAMPLES 

load 

3  1 

1  8        10.         0. 
load 

1  1 

1  6       -10.  2 


APPENDIX  D 

P5-NFAP 

FE  PROGRAM  TO  ANALYZE  PLANE  STRAIN  PROBLEMS 

WITH  THE  CAP  SOIL  MODEL 
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INTRODUCTION 

PS-NFAP  Is  a  finite  element  (FE)  computer  program  for 
analysis  of  plane  strain  problems  using  the  cap  strain- 
hardening  model  to  represent  soil  behavior.   It  Is  capable 
of  analyzing  embankments  constructed  on  soft  ground  and  has 
•pec la  1  features  which  allow  analysis  of  reinforced  ambank- 
aients.   PS-NFAP  was  adapted  by  Humphrey  (1965)  from  a  gen- 
eral purpose  FE  program  called  NFAP  that  was  originally  de- 
veloped by  Chang  (I960).   The  cap  model  was  Implemented  In 
NFAP  by  HcCarron  (19B5)  and  McCarron  and  Chen  (1986a. 
1986b).   PS-NFAP  Is  written  primarily  In  FORTRAN  IV  although 
some  sections  of  the  code  use  features  of  FORTRAN77.   It  Is 
Implemented  on  an  IBM-PC  using  Ryan- McFar land  FORTRAN,  ver- 
sion 2.0  ( Ryan- McFar land,  1985)  and  on  an  IBM  3083  mainframe 
using  VS-FORTRAN.   The  program  uses  the  convention  that  com- 
pressive stresses  and  strains  are  negative.   Since  most 
stresses  encountered  In  soli  mechanics  are  compressive,  most 
stresses  for  analysts  of  soils  problems  with  PS-NFAP  will  be 
negative. 

This  appendix  is  organized  as  follows.   The  general 
capabilities  of  PS-NFAP  are   described  In  the  next  section; 
then  program  organization  Is  outlined.   Next  program  opera- 
tion Is  described.   Input  instructions  are   given  In  Appendix 
0.1  and  comments  to  the  input  instructions  are  given  in 
Appendix  0.2.   An  example  illustrating  use  of  PS-NFAP  to 
analyze  a  reinforced  embankment  constructed  on  soft  ground 
Is  given  In  Appendix  E. 

CAPABILITIES  OF  PROGRAM 

PS-NFAP  performs  an  incremental  load-displacement  anal- 
ysis.  After  each  increment  of  load  Is  applied  the  displace- 
ment field  Is  modified  ustng  an  Iterative  procedure  until  an 
equilibrium  configuration  Is  reached.   Convergence  Is  based 
on  the  difference  between  two  successive  displacement  norms 
as  defined  in  Section  5.2  of  Humphrey  and  Holtz  (1986). 

Two  element  types  are    Included  In  PS-NFAP.   The  first 
Is  a  two  or  three  node  TRUSS  element  which  can  carry  only 
axial  load.   Material  behavior  is  either  linear  elastic  or 
nonlinear  elastic.   This  element  is  used  primarily  to  repre- 
sent tensile  reinforcement.   There  is  an  option  that  allows 
the  element  to  be  initially  Inactive  and  then  subsequently 
activated  at  a  specified  load  step.   This  can  be  used  to 
simulate  one  or  more  reinforcing  layers  placed  In  an  embank- 
ment during  construction  at  a  level  above  the  original 
ground  surface. 

The  second  element  type  Is  a  two-dimensional  (2-0)*  4 
to  8-node,  Isoparametric,  plane  strain  continuum  element. 
It  Is  used  mainly  to  represent  the  foundation  soil  and 
embankment  fill.   Soil  behavior  Is  represented  by  the  cap 
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model.   The  problem  can  be  formulated  considering  only 
material  nonllnearity  or  with  an  updated  Lagrenglen  formula- 
tion.  The  latter  expresses  equilibrium  equations  fn  terms 
of  the  most  recent  equilibrium  configuration  and  accounts 
for  both  geometric  (large  displacement)  and  material  non lin- 
earities (HcCarron,  19BS).   Details  of  the  cap  model  are 
given  In  HcCarron  (1985)  and  Chen  and  Baladl  (1985)  and  pro- 
cedures to  determine  the  cap  parameters  are   given  fn 
Humphrey  and  Holtz  (1986).   The  model  allows  Initial 
stresses  to  be  specified.   The  element  can  be  Initially 
Inactive  and  then  activated  at  a  specified  load  step.   This 
Is  useful  to  simulate  lift  by  lift  embankment  construction. 
Behavior  can  be  drained  or  undralned.   There  Is  an  option 
which  allows  behavior  to  be  switched  from  drained  to 
undralned  at  a  specified  load  step.   This  has  been  used  to 
establish  the  Initial  state  of  stress  for  an  existing 
embankment  on  a  fully  consolidated  foundation  prior  to 
placement  of  additional  fill  under  undralned  conditions 
(Humphrey  and  Holtz,  1986).   The  2D  element  also  has  a 
linear  elastic  material  behavior  model. 

In  addition.  PS-NFAP  has  a  routine  to  renumber  the 
nodes  to  minimize  the  band  width  of  the  stiffness  matrix. 
thereby.  Improving  solution  efficiency.   In  the  IBM-PC  vei — 
sfon,  this  is  Implemented  as  a  separate  program  called 
RENUMX.   There  Is  also  an  out  of  core  equation  solver  that 
permits  solution  of  large  problems  on  computers  with  limited 
memory.   Solution  results  can  be  saved  for  use  as  the  start- 
ing point  of  a  restart  Job  in  which  additional  load  steps 
are  appl ted. 

PROGRAM  ORGANIZATION 

PS-NFAP  is  comprised  of  a  MAIN  program  that  calls 
subroutine  NFAP  which  controls  the  remainder  of  the 
execution.   The  57  subroutines  used  In  the  main  portion  of 
PS-NFAP  ere  listed  In  Table  D.l  along  with  a  short  descrip- 
tion of  their  purpose.   The  subroutines  used  in  the  node 
renumbering  portion  of  the  program  are    listed  In  Table  D.2. 
There  Is  also  a  BLOCK  DATA  unit. 

Program  operation  can  be  separated  into  input,  solu- 
tion, and  equilibrium  Iteration  phases.   An  outline  of  the 
operations  performed  In  each  phase  Is  given  in  HcCarron  and 
Chen  (1986a).   A  simplified  flow  chart  of  PS-NFAP  operation 
Is  shown  In  Fig.  D.l.   Some  of  the  Input  data.  Intermediate 
calculations,  and  solution  results  are   stored  in  a  large 
matrix  called  the  *A-matrlx'  which  Is  assigned  to  blank 
common.   Other  data  Is  temporarily  stored  In  files  for  use 
during  subsequent  load  steps.   File  usage  Is  summarized  In 
Table  D.3. 
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PROGRAM  OPERATION 

Input  Instructions  are  given  In  Appendix  D.I  and  com- 
ments to  the  Input  Instructions  are  given  In  Appendix  0.2. 
The  Input  file  Is  assigned  to  standard  Input  (unit  5)  and 
results  mre   output  to  standard  output  (unit  6).   Example 
Input  and  output  files  are  shown  In  Appendices  E.3  and  E.5. 
respectively. 

IBM-PC  Version 

The  IBM-PC  version  runs  on  IBM-PC,  IBM-XT,  IBM-AT,  or 
compatible  computers  with  at  least  512K  of  memory  and  a  hard 
disk.   More  memory,  up  to  the  system  limit  of  640K,  Is 
desirable  for  larger  problems.   The  first  step  In  running 
PS-NFAP,  If  the  node  renumbering  option  is  chosen.  Is  to  run 
the  node  renumbering  program  NFMINX.   It  Is  executed  with 
the  fol lowing  command: 

ONFMINX  /R  20000  <F1LE.NFI  >FILE.NM0 

where  C>  is  the  PC-DOS  prompt,  /R  20000  reserves  20K  of  mem- 
ory for  records  from  unformatted  read  and  write  statements 
(Ryan-McFarland,  1985),  FILE.NFI  is  the  PS-NFAP  Input  file. 
and  FILE.NMO  Is  the  output  file.   A  list  of  the  renumbered 
nodes  is  written  to  FILE.NMO  and  to  a  file  named  NFMIN.NFP 
for  later  use  by  PS-NFAP.   An  example  output  Is  shown  In 
Appendix  E.4. 

Next  PS-NFAP  is  Invoked  with  the  following  command: 

OPSNFAP  /R  63735  <F1LE.NFI  >FILE.NF0 

where  /R  63735  reserves  63. 7K  of  memory  for  records  from 
unformatted  read  and  write  statements  (Ryan-McFarland,  1985) 
and  FILE.NFO  1s  the  output  file.   An  example  output  Is  shown 
in  Appendix  E.5.   For  large  problems  the  temporary  scratch 
files  listed  In  Table  D.3  and  FILE.NFO  require  considerable 
disk  storage  In  some  cases  exceeding  1  megabyte.   The  user 
should  check  that  there  Is  sufficient  disk  storage  prior  to 
execution  of  PS-NFAP. 

If  It  Is  necessary  to  recompile  PS-NFAP  using  the  Ryan- 
HcFarland  (1985)  compiler  the  /i  option  should  not  be  used. 
This  causes  the  default  size  of  all  integers  to  be  INTEGER*4 
(32  bits).   Subroutine  ASSEM  only  should  be  compiled  with 
the  /b  option  to  generate  code  to  access  arrays  larger  than 
64K. 

Size  of  A-matrlx 

Large  problems  can  be  run  most  efficiently  If  there  Is 
sufficient  space  In  the  A-matrtx  to  allow  the  solution  to  be 
performed  completely  In  core.   However,  the  maximum  size  of 
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the  A-matrlx  may  be  limited  by  the  computer  system.   On  an 
IBM-PC  with  640K  of  memory  the  maximum  size  Is  72000. 

The  size  of  the  A-matrlx  can  be  changed  If  necessary. 
In  the  mainframe  version  this  Is  done  by  changing  the  size 
of  the  A-matrlx  and  the  variable  MTOT  in  the  main  program, 
recompiling  the  main  program,  and  then  relinking  the  pro- 
gram.  Note  that  the  size  of  the  A-matrlx  Is  MTOT  4  10.   In 
the  IBM-PC  version  the  size  Is  changed  by  altering  the  size 
of  A,  I A,  and  MTOT  in  the  file  named  A.CMM  which  is  shown 
below: 

COMMON  A(72010) 
DIMENSION  1A(72010) 
EQUIVALENCE  (A(l).IA(l)) 
MTOT  «  72000 

Note  that  the  size  of  A  and  1A  are  10  greater  than  the  size 
of  MTOT.   The  following  program  units  which  contain  the  A- 
metrlx  are   then  recompiled:  EL2D10,  ELCAL.  LSTM,  MAIN,  NFAP, 
NFAPIN,  RUSS,  STRCAL,  TDFE,  TODMFE .  and  TRUSS  using  the 
batch  file  RECOMA.BAT  listed  In  Table  D.4.   The  entire  pro- 
gram Is  then  linked  using  the  link  file  PSNFAP.LNK  listed  In 
Table  D.5  in  conjunction  with  the  Ryan-McFar land  linker 
using  the  command  shown  below  (Ryan-McFar land,  1985): 

OPLINK86  GPSNFAP 
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Tab 1 e  D . I 
PS-NFAP  subroutines 


ACTRED     Equation  solution 

ADDBan    Save  element  stiffness  on  tape  &  assemble  load 

vector 
ADDMAT    Assemble  element  stiffness  matrices  Into  Qlobal 

stiffness  matrix 
ADDRES    Calculate  addresses  of  diagonal  elements  In  banded 

matrix 
ASBLK      Assemble  element  blocks 
ASOLVE     Equation  solution;  Cnoleskl  decomposition  of 

profile  stored  In  symmetric  matrix 
ASSEM     Assembles  nonlinear  stiffness  matrix 
BACSUB     Equation  solution 

CAUCHY     Convert  Plola-KI rkhhof  stresses  to  Ceuchy  stresses 
COLHT     Computes  column  height  of  banded  matrix 
DERIQ     Evaluates  strain-displacement  matrix  for 

quadrilateral  element 
DEVPRT     Computes  devletorlc  stresses 
D0T1       Computes  vector  dot  product 
EL2D10     CAP  model 
ELCAL      Input  element  data 
ELEMNT     Calls  routines  for  reading,  generating  and  storing 

element  data 
EP2D10     CAP  model 
EPL2  10     CAP  model 

EQUIT      Equation  solution;  equilibrium  Iteration 
FUNCT2     Finds  Interpolation  functions  and  derivatives  for 

4-  to  8-node  Isoparametric  quadrilateral;  finds 

Jacoblan  and  Its  determinant 
FUNCTC     Finds  Interpolation  functions 

truss  elements 
GRAV       Computes  gravity  force  during 

embankment  construction 
1EP210     Initializes  working  array   for 
IN1TAL     Saves  ID  array   and  Initializes 

zero 
INITWA    Calls  subroutines  to  Initialize  working  array  for 

2  dimensional  material  models 
INPUT      Input,  generate  and  print  nodal  data;  compute 

equation  numbers  (ID  array) 
LOADEF     Calculate  effective  loads  for  nonlinear  elements 

and  updates  load  vector  calculated  In  LOADVC 
LOADS     Reads  load  curves  and  nodal  load  data;  calculates 

load  vector  for  each  time  step 
LOADVC    Calls  NFAPIN  to  read  and  calculate  load  vectors 
LSOLVE     Equation  solution 
LSTM      Assembles  linear  stiffness 


for  2-  or  3- node 
1 ncremental 


CAP  model 
dl splecements 


to 


iiittiiimi 
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PS-NFAP 


Table  D. J 
subroutines 


( cont . ) 


BBSSBEEEEeEEEBBEEBEEBBBEKBBEBBSBBEZEEEEEEBEBBBEBBBBBKBKSEEEZ 

HATBAR    Material  model  for  truss  element 

MATRT2    Print  materiel  properties  for  quadrilateral 

element 
MAXMJN    Computes  principal  stresses 
MD2D10    CAP  model;  develops  elastic-plastic  stress-strain 

relation  If  loading  occurs  on  CAP 
HKBLK     Sets  up  blocks  for  equation  storage 
MULT      Equation  solution;  computes  product  of  blocked 

matrix  and  vector 
NEWDAV    Calculates  new  displacements  at  time  T+delta(T) 
NFAPIN     Reads  Input  data 
PF1LE      Equation  solution;  reads  and  writes  records  to 

tape 
PLACE      Updates  current  displacement 
PRAGER    Forms  elasto-plest 1 c  material  matrix  for  Drucker- 

Prager    yield  criterion 
QUADC     Computes  element  stiffness  and  stresses  for  truss 

element 
QUADS     Computes  element  stiffness  for  Isoparametric 

quadr I  1 atera 1  e 1 ement 
RSTART    Saves  data  for  restart  Job 
RUSS       Truss  element 

SECOND     Monitors  elapsed  time  (IBM-PC  version  only) 
SIZE      Check  of  storage  requirements 
STRCAL     Calculate  maximum  block  size  and  determine  If  In 

core  or  out-of-core  solution  scheme  must  be  used 
STRESS    Calculates  nodal  stresses 
STSTL     Generate  global  stress-strain  law  for  Isotropic, 

linear  materials  In  plane  strain 
STSTN     Generate  global  stress-strain  law  and  stresses  for 

nonlinear  material  models 
SUMLD     Computes  and  prints  loading  function  table 
TDFE      Quadrilateral  element 
TODMFE     Quadrilateral  element 
TRUSS     Truss  element 
WRITE      Prints  displacements 


-  DB 


Table  0.2 
Subroutines  used  for 
nodal  numbering  optimization 


BAND 

CHECK 

DELETE 

DGREE 

ELSTR 

FNDIAM 

FORMLV 

G1BB 

LINK 

NFPMJN 


NINCON 

NMBER 

P1KLVL 

PROF  1 T 

REDUCE 

SETAP 

SORT2 

SORTDG 

TREE 


Table  D.3 
File  usage  by  PS-NFAP 


File   Subroutines 
no.   referenced 


Information  stored  fn 
f  1  le 


1  ASSEM,ELCAL, STRESS 

2  ASSEM,ELCAL,EQU1T, 
RSTART .STRESS 

3  EQUIT.LOADEF, LOADS. 
NFAP 

4  ADDBAN.ASBLK.E0U1T. 
MULT, NFAP, PFILE 


8   INI TAL. LOADS. NFAP. 
RSTART, WRITE 


9  ASSEM.EQU1T, RSTART, 
STRESS 

10  EQU1T, NFAP, PFILE 

12  ADDBAN.ASBLK, PFILE 

14  MKBLK, PFILE 

18  PFILE 

20  NFAPIN.GIBB 


Note:  Files  1,  2,  3,  4,  10 

and  are  not  saved  at 


Linear  element  group  data 
Nonlinear  element  group  data 

Externally  applied  loads 

(I.)  Linear  element  stiffness 
matr lx, element  by  element 

(2.)  Linear  structural  stiffness 
in  blocks 

(1.)  ID  array  and  Initial  dis- 
placement 

(2.)  ID  array,  displacement  and 
displacement  Increments  for 
restart 

File  Is  named  TAPES 

Same  as  file  2  and  nonlinear 

group  data  for  restart 

File  Is  named  TAPE9 

Effective  nonlinear  stiffness  In 

blocks 

Global  stiffness  matrix  in  blocks 

Designated  block  structure  for 

assemblage 

Element  stiffness  before 

assemblage 

Renumbered  nodes  In  file  named 

NFMIN.NFP  (IBM-PC  ver»lon  only) 

IIIItltllKIttltlltUUMRUIttlt* 

,  12,  14,  and  18  are   scratch  files 
the  end  of  program  execution. 
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Table  D.4 
Batch  file  RECOMA.BAT  to  recompile  subroutines 
when  size  of  A-matrix  Is  modified 


:sssb  =  s  =  sb  =  : 


REM    BATCH  FILE  TO  RECOMPILE  SUBROUTINES  WHEN  SIZE  OF 

REM    A- MATRIX  MODIFIED 

REM 

RMFORT  EL2D10  /L  >EL2D10.LST 

RMFORT  ELCAL   /L  >ELCAL.LST 

RMFORT  LSTM    /L  >LSTM.LST 

RMFORT  PSNFAP  /L  >PSNFAP.LST 

RMFORT  NFAPS   /L  >NFAPS.LST 

RMFORT  NFAPIN  /L  >NFAPIN.LST 

RMFORT  RUSS    /L  >RUSS.LST 

RMFORT  STRCAL  /L  >STRCAL.LST 

RMFORT  TDFE    /L  >TDFE.LST 

RMFORT  TODMFE  /L  >TODMFE.LST 

RMFORT  TRUSS   /L  >TRUSS.LST 


Table  D.5 
Link  file  PSNFAP. LNK  to  link  PS-NFAP 


F I LE  PSNFAP , BLKDATA , NFAPS , RSTART .SIZE, SECOND 

FILE  ADDRES, ELCAL, ELEMNT, JNITAL, INPUT, LOADS, LOAD VC 

FILE  LSTM , MKBLK , NF AP I N , STRCAL , ACTRED 

FILE  ADDBAN . ADDMAT , ASBLK , ASOL VE , ASSEM , BACSUB . COLHT , DER I Q 

F I LE  DOT  I , EQU I T , LOADEF , LSOLVE , MULT , NEWDAV , PF I LE , PLACE , STRESS 

FILE  SUMLD, WRITE 

F I LE  CAUCHY , I N I  TWA , MATRT2 , MAXM 1 N , QUADS , STSTL , STSTN 

FILE  TDFE. TODMFE 

F I LE  DE VPRT , EL2D 1 0 , EP2D 1 0 , EPL2 1 0 . FUNCT2 . GRA V , I EP2 I  0 

FILE  MD2D10.PRAGER 

FILE  FUNCTC . MATBAR .QUADC , RUSS , TRUSS 

OUTPUT  PSNFAP 
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Figure  D.I.   PS-NFAP  flowchart. 


APPENDIX  D.l 
INPUT  INSTRUCTION 


-  Dll  - 


EAD1NG 
ard  1.1 

oluint    Variable 
-  72     HED  (18) 


CoBnent*  Entry 

(1)         Enter  the  master  heading  information 
for  use  In  labeling  the  output 


Detailed  comment  on  the  definition  of  variable  can  be  found  in 
APPENDIX  D.2. 
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MASTER  CONTROL  CARDS 

Card  2.1 

Columns    Variable    Comment 

1-5        NUMNP        (1) 


12-15 


16-20 


21-25 


36-40 


A6-50 


51-55 


NEGL 


NEGNL 


MODEX 


26-30     •  NPER 
31-35      1PR1 


IRINT 


1STOTE 


NRENUM 


(2) 


(3) 


(A) 


(5) 
(6) 

(7) 
(6) 
(9) 


Entry 

Total  nutter  of  nodal  point*; 
EQ.O;  program  stop 

Number  of  linear  element  groups; 
EQ.O;  all  elements  are  nonlinear 

Number  of  nonlinear  element  groups; 
EQ.O;  all  elements  are  linear 

Flag  Indicating  solution  mode; 

EQ.O 

EQ.l 

EQ.2 


data  check 

execution 

restart 


Number  of  different  solution  time  periods 

Output  printing  Interval; 
EQ.O;  default  set  to  "1" 

Restart  save  time  step  Interval 
EQ.O;  no  saving  for  restart 
GT.O;  Interval  for  data  saving 

Specified  block  size 

EQ.O;  calculated  automatically 

by  the  program 

Flag  Indicating  node  renumbering  to 

minimize  column  height  of  stiffness 

matrix 

EQ.O;  no  node  renumbering 

EQ.l;  node  renumbering 
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MASTER  CONTROL  CARDS  (continued) 
Card  2.2 

Columns   Variable   Comment  Entry 

1-5        1SREF      (1)       Number  of  steps  between  re-for»lng 

effective  itiffness  matrix 
EQ.O;  default  set  to  1 

6-10       NUMREF     (2)       Number  of  Iterations  between 

stiffness  reformation  each  time 

a  tep 

EQ.O;  no  reformation 

11-15      1EQU1T     (3)       Number  of  time  steps  between  equili- 
brium iterations 
EQ.O;  default  set  to  1 

16-20      ITEMAX     (3)       Maximum  number  of  equilibrium  Itera- 
tions permitted  in  each  time  step 
EQ.O;  default  set  to  15 

21-25      IACC  Flag  indicated  on  acceleration  scheme 

used  during  equilibrium  iteration 
EQ.O;  no  acceleration 
EQ.l;  secant  acceleration 

41-45      RTOL       (3)       Relative  tolerance  used  to  measure 

displacement  convergence 
EQ.O;  default  set  to  "l.E-3" 

46-50     FTOL        (3)       Relative  force  tolerance  used  to 

measure  equilibrium  convergence 
EQ.O;  default  set  to  "0". 
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MASTER  CONTROL  CARDS  (continued) 

Card  2.3 

CoIubds  Variable        Comment 

1-10  TSPER(l)          (1) 

11-20  DTPER(l)          (1) 

21-30  TSPER(2) 

31-^0  DTPER(2) 


Entry 
Starting  tine  for  1st  period 
Tine  atep  increment  for  lat  period 
Starting  tine  for  2nd  period 
Tine  atep  increment  for  2nd  period 


71-80 
Next  card 

1-10 
1-20 


DTPER(4) 

(If  required) 

TSPER(5) 

DTPER(5) 


Tine  atep  Increment  for  4th  period 

Starting  tine  for  5th  period 

Tine  atep  Increment  for  5th  period 


TFINAL 


(2) 


Final    time    of    analysis 
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AkSTEk  CONTROL  CARDS  (continued) 
Card  2.4 

ftoluint    Variable 
RPB 


P-10 


1-25 


IDC 


IDSET 


Comment 
(1) 


Ent  ry 


(2) 


(3) 


Number  of  blocks  of  displacement 

print  out. 

EQ.O;  print  all  nodal  point 

c  omponents 

Displacement  printout  code 
EQ.O;  no  displacement  printout 
EQ.l;  print  displacements 

A  reference  node  for  displacement 
control  of  ghost  elements 
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MASTER  CONTROL  CARDS  (continued) 
Skip  this  card  If  NPB.EQ.O 

Card  2.5 

Variable      Comment 
IPNODE(l.l)    (1) 
IPNODE(2,l ) 


ColUBDS 

1-5 

6-10 

11-15 


1PNODEU.2) 
etc. 


Ent  ry 
First  node  of  printout  block  No.  1 
Last  node  of  printout  block  No.  1 
First  node  of  printout  block  No.  2 
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DAL  POINT  DATA 

(Input  as  many  cards  as  required  for  all  nodes) 
Holuans    Variable    Comment 


L-5 


.-70 


Entry 


K 


PSF 


(1) 
(2) 


ID(2,K) 
IDO.K) 

Y(N) 
Z(N) 

KN 


(3) 


(4) 


Node  (joint)  number: 
CE.l  and  LE.NUMNP 

Print  suppression  flag  (ignored  unless 

N.EQ.l); 

EQ.  ;  (blank.)  no  suppression 

EQ.A;  suppress  the  list  of  nodal  coordinates 

EQ.B;  suppress  the  list  of  equation  numbers 

EQ.C;  both  A  and  B  above 

EQ.D;  suppress  the  list  of  generated  nodal 

coordinates 

Y-t r ans 1  at  ion  boundary  code 
Z-t ranslat Ion  boundary  code 

Y-coordinate 
Z-coordinat e 

Node  number  Increment  for  node  data 

genera t  ion ; 

EQ.O;  no  generation 


-  Die  - 


LOAD  CONTROL  CARDS 
Card  4.1 

Columns    Variable 
1-5        KLOAD 


6-10 


11-15 


ULCUR 


KPTM 


Connent 
(1) 

(2) 

(2) 


Entry 

Number  of  nodal  force  components  to 
apeclfy  nodal  forcea 

Total  number  of  load  curves 
( t lne  f unct Ions ) 

Maximum  nunber  of  polnta  used  to 
describe  any  one  of  the  load  curves 


41-45 


NTURNS 


(3) 


Nunber  of  reverae  load 


1 
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OAD  CONTROL  CARDS  (continued) 


Skip  this  card  If  NTURNS.EQ.O 
f.rd  4.2 
Coluant    Variable    Consent 


f 


-10 

■  1-20 

21-30 


TRV( 1) 
TRV(2) 
TRV(3) 


(1) 


Ent  ry 
Tine  of  the  first  load  reverae  point 
Tine  of  the  aecond  load  reverae  point 
Tine  of  the  aecond  load  reverae  point 


i 

I 
I 

■ 
I 

1 

1 

I 
1 
1 
I 
I 


1-80      TRV(6)  Time  of  the  eighth  load  reverse  point 

oput  NTURNS  points;  use  additional  cards  as  required. 


-  D20  - 


ELEMENT  LIBRARY 

Element  type  Is  Identified  by  the  first  entry  (NPAR(l))  of 
the  first  card  In  this  section,  i.e., 


Element 


NPAR(1  ) 


Truss  1 

2-D  element      2 


A  structure  say  be  divided  Into  several  groups  of  elements, 
consisting  of  either  linear  groups,  nonlinear  groups  or  both. 
Input  as  sany  blocks  of  data  In  this  section  ss  there  are  total 
eleaeot  groups.  Total  number  of  linear  element  groups  (NEGL)  are 
Input  first,  and  followed  by  the  total  number  of  nonlinear  ele- 
ment groups  (NEGNL). 

In  any  one  group  all  element  Input  must  be  the  same  element 
type  and  analysis  type;  e.g.  If  nonlinear  TRUSS  elements  are 
given  as  input,  then  all  elements  In  this  group  must  be  non- 
linear.  Furthermore,  In  any  one  group,  only  one  material  model 
can  be  used,  e.g.,  all  elements  In  the  group  must  be  defined  by 
the  cap  model.   However,  a  number  of  different  sets  of  material 
constants  for  a  specific  model  can  be  used. 
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^RUSS  ELEMENT 

«RUSS  tleientt  are  two  or  three-node  ■embers  allowed  arbitrary  orientation 
n  the  Y,  Z  aystem.   The  TRUSS  transmits  axial  force  only,  and 
In  general  Is  a  alx  (6)  degree  of  freedom  element  (I.e.,  two 
jalobal  translation  components  at  each  of  the  element  nodes). 


Card    5.1 

foluint 
1-4 

Is 

f 

T-12 


6-60 


il-64 


5-68 


Variable 
NPAR(  1  ) 
NPAR(2) 

NPAR(3) 


NPAR(15) 


NPAR(16) 


Cosiment 


(1) 


(2) 


(3) 


(A) 


NPAR(17)    (5) 


Entry 

Enter  the  number  "1" 

Number  of  TRUSS  elements  In 

this  group : 

CE.l 

Type  of  analysis 

EQ.O;  linear  analysis 

EQ.l;  materially  nonlinear  only 

EQ.2;  large  displacement 

Material  model  number 
EQ.l ;  linear  elastic 
EQ.2;  nonlinear  elastic 

Number  of  different  sets  of 
section/material  properties 
GE.l 

Number  of  material  model 
constants  per  aet 
EQ.l  If  NPAR(15).EQ.l 
GE.4  If  NPAR(15).EQ.2 
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TRUSS  ELEMENT  (continued) 

Input  cards  for  linear  elastic  aaterlal  and  truss  member  properties. 
Two  loput  cards  are  required.   Skip  this  aet  of  cards  If  HPAR( 1 5 ) . NE . 1 
Otherwise  read  NPAR(16)  aets  of  cards. 


Card  5.2a.   Material  number  card 
Columns    Variable    Comment 
1-5        K 


En  t  ry 


Material /sec t Ion  number 
GE.l ;  and  LE.NPAR( 16) 


Card  5.2b.   Property  card 
Columns    Variable    Comment 
1-10       E(N)        (1) 
11-20      AREA(N) 
31-40      PINIT(N) 


Ent  ry 
Young's  nodulus 
CroEE-sectlonal  area 
Initial  axial  force 


1 

I 
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RUSS  ELEMENT  (continued) 

tput  cards  for  nonlinear  elastic  material  and  ■ember  property 
rds.   Two  or  tore  Input  cards  are  required.   Skip  this  aet  of 
carda  If  NPAR( 1 5) .NE.2 ,  otherwise  read  NPAR(16)  lets  of  cards. 


rd  5.3a.   Material  number  card 


rr 


Variable 


Comment 


K 


t 


r  1-20 
■1-60 
41-50 


AREA(N) 

PINIT(N) 

TIME-ON 


(1) 
(2) 
(3) 


En  t  ry 

Mat erlal /aec t Ion  number 
CE.l  and  LE.NPAR(16) 

Cross-sectional  area 

Initial  axial  force 

Time  element  is  activated 


I 


Card  5.3b.   St ress-st rain  curve  card  (8F10.0) 
loluonB        Variable        Comment 


■-10 
11-20 


Ent  ry 


PR0P(  1  ,K) 
PR0P(2,N) 


(4) 


Strain  at  point  1,  c 
Strain  at  point  2,  c' 


1 


PR0P(NC0N/2,N) 

PR0P(NC0N/2+l,N) 

PR0P(NC0N/2+2,N) 


Strain  at  point  NCON/2 

Stress  at  point  1,  o 

2 
Stress  at  point  2,  o 


PROP(NCON.N) 


Stress  at  point  NCON/2 
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TRUSS  ELEMENT  (continued) 

Card  5.4 

MPAR(2)  elements  Bust  be  Input  and/or  generated  In  this  section 
In  ascending  sequence  beginning  with  "1". 


Columns 
1-5 

6-10 
11-15 
16-20 
21-25 

26-30 


31-35 


Variable 
M 

II 
JJ 
KK 

MTYP 

IPS 


Comment 


KG 


(1) 
(1) 
(1) 


(2) 


Ent  ry 

TRUSS  element  number; 
CE.l  and  LE.NPAR(2) 

First  Node 

Second  Node 

Third  Node 

Material  property  set  number; 
GE.l  and  LE.NPAR(16) 

Flag  for  printing  axial  stress 

in  TRUSS  element; 

EQ. 1 ,  no  printing 

EQ.O,  print  element  stress 

Node  generation  Increment  used 
to  compute  node  numbers  for 
missing  elements; 
EQ.O;  default  set  to  "1" 
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J-D  CONTINUUM  ELEMENT 

trd  6.1 
Iuidi    Variable    Comment 
f-4        WPAR(l) 

HPAR(2)     (1) 


L-12 


.-28 


'-40 


r-60 


-64 


.-68 


i-72 


NPAR(3)     (2) 


NPAR(7)     (3) 


NPAR(IO)    (4) 


NPAR(  15) 


NPAR(16) 


(5) 


(5) 


NPAR(17)    (5) 


NPAR(18)    (5) 


Entry 

Enter  the  number  "2" 

Number  of  2-D  CONTINUUM  elements 

In  t h  i  6  group: 

CE.l 

Flag  indicating  type  of  analysis 

EQ.O;  linear  analysis 

EQ.l;  material  nonlinear  analysis  only 

EQ.4;  updated  Lagrangian  (large  rotation) 

Maximum  number  of  nodes  used  to 
describe  any  one  element: 
CE.4  and  LE.8 
EQ.O;  default  set  to  "8" 

Numerical  Integration  order  to  be 
used  In  Gauss  quadrature  formula: 
EQ.O;  default  set  to  "2" 
GE.2  and  LE.4 

Material  model  numbers: 
EQ.l;  linear  isotropic 

EQ.10;  cap  model  with  Dr ucker-Prager  Shear 
limitation 

Number  of  different  sets  of  material 

properties : 

CE.l 

Number  of  constants  per  property  set; 
EQ.O;  if  NPAR(15).EQ.10 
EQ.22;if  NPAR(15).EQ.10 

Dimension  of  storage  array  required 
for  element  history; 
EQ.O;  if  NPAR(15).EQ.l 
EQ.12;if  NPAR(15).EQ.10 
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2-D-CONTINUUM  ELEMENT 

NPAR()6)  sets  of  Card  6.2  and  Card  6.3  Bust  be  input.   Card  6.2  1b  the 
sane  for  all  aaterlal  aodels,  but  Card  6.3  depends  on  the  uteris! 
■odel  Duaber  (NPAR(15));  choose  from  Card  6.3.1  or  Card  6.3.10. 

Card  6.2 


Columns    Variable    Comment 
1-5        N 


Entry 

Material  property  aet  number; 
GE.l  and  LE.NPAR(16) 
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2-D  CONTINUUM  ELEMENT 

Bard  6.3.1  -  Linear  isotropic  -  NPAR( ] 5 ) . EQ . 1 

toluana    Variable     Comment  Entry 

^1-10      PROP(l.N)    (1)  Young's  »odulus,  E 

■11-20     PROP(2,N)  Polsson's  ratio,  v 
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2-D  CONTINUUM  ELEMENT 

Card  6.3.10  -  Soil  Plasticity  Cap  Model 
NPAR(15).EQ.10 

Each  Card  Contains  8  Input  data.  (8F10.0) 


Columns 
Card  8.3.10a 
1  -  10 
1  1  -  20 

21  -  30 
31  -  40 
41-50 
51  -  60 
61  -  70 
71-80 


Variable    Comment    Entry 


PR0P(  1) 
PR0P(2) 
PROP(3) 
PR0P(4) 
PR0P(5) 
PR0P(6) 
PR0P(  7) 
PR0P(8) 


(1 
(1 
(1 
(1 
(1 
(1 
(2 
(2 


nln 


1 


2 
AM 

AC 


Card  6.3.10b 

1  -  10 

PR0P(9) 

(3) 

AW 

1  1  -  20 

PR0P(  10) 

(3) 

AD 

21-30 

PR0P( 1 1) 

(4) 

R 

31  -  40 

PR0P(  12) 

(5) 

XL 

41  -  50 

PR0P(13) 

(6) 

TENCUT 

51  -  60 

PR0P(14) 

(7) 

STATE 

61-70 

PR0P(15) 

(8) 

Al 

71  -  80 

PR0P(16) 

(8) 

A, 

I 
I 

t 

ColuBOB 

fard  6.3.10c 

i,"0 

VI  -  20 

■l  -  30 

31  -  AO 

§  -  50 

1  -  60 
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t 


D  CONTINUUM  ELEMENT 
rd  6.3.10   (Continued) 
Variable 

PR0P(17) 
PR0P(18) 
PR0P(19) 
PROP(20) 
PROP(21) 
PROP(22) 


Comment 

(8) 
(9) 

(10) 
(11) 
(12) 
(12) 


Ent  ry 

K0 
SWITCH 

PAC 

TIME-ON 

WCT 

BOUY 
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2-D  CONTINUUM  ELEMENT 
Card  6.4a 

(Two  cards  are  needed  for  each  element.) 
Columns    Variable    Consent 
M 


1  -  5 


6-10 


16-20 


21  -  25 


Entry 


IEL 


MTYP 


KG 


(1) 
(2) 


(1) 


2-D  CONTINUUM  element  nunber; 
CE.l  and  LE.NPAR(2) 

Nunber  of  nodes  used  to  describe 

this  elenent ; 

EQ.O;  default  set  to  "NPAR(7)" 

LE.NPAR(7) 

Material  property  set  nunber 
assigned  to  this  elenent; 
GE.l  and  LE.NPAR(16) 

Node  generation  parameter  used  to 
conpute  node  nunbers  for  missing 
elements  (given  on  first  card  of 
a  sequence ) ; 
EQ.O;  default  set  to  "1" 
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D  CONTINUUM  ELEMENT  (continued) 
rd  6.4b 

Iuidi    Variable    Comment 
W0D(1)      (2) 


-  10 


I 


r 
r 


Ent  ry 


K0D(2) 
i  -  15    N0D(3) 


16  -  20    N0D(4) 

I 


21  -  25    N0D(5) 
6-30    N0D(6) 


1  -  35    NOD(7) 


-  40    N0D(8) 


Global  node  number  of  element 

nodal  point  1 

Global  node  number  of  element 

nodal  point  2 

Global  node  number  of  element 

nodal  point  3 

Global  node  number  of  element 

nodal  point  4 

Global  node  number  of  element 

nodal  point  5 

Global  node  number  of  element 

nodal  point  6 

Global  node  number  of  element 

nodal  point  7 

Global  node  number  of  element 

nodal  point  8 
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APPLIED  LOADS 

Supply  NLCl'R  lets  with  2  or  sore  cards  per  set. 

Card  7.U  (215) 


Columns    Variable    Comment 
1-5      NTF 


6-10 


HPTS 


En  t  ry 

Tine  function  nuaber; 
CE.l  and  LE.NLCUR 

Number  of  points  (i.e.,  f(t), 

t  pairs)  used  to  input  this  time 

function; 

CE.2  and  LE.NPTM 


Card  7.1b 

ColunnE   '  Variable 


1  -  10 

1  1  -  20 

2  1  -  30 
31  -  40 


TIMV( 1) 
RV(1) 
TIMV(2) 
RV(2) 


Comment  Entry 

(1)        Time  at  point  1,  t 


1 


Function  at  point  1,  f(t.) 

Time  at  point  2,  t- 

Functlon  value  at  point  2,  f(t  ) 


71  -  80  RV(4) 

Next  card  (if  required) 

I  -  10  TXMV(5)    (2) 

II  -  20  RV(5) 


Function  value  at  point  4,  f(t.) 


Time  at  point  5,  t 

Function  value  at  point  5,  f(t  ) 


I 

! 
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PPLIED  LOADS 
■kip  this  section  of  NLOAD.EQ.O 
Card  7.2 

j^olumnt    Variable    Comment 
-  5      MOD         (1) 


-  10 


I 
| 
I 

11-15 

I 


16  -  25 


1DIRN 


NCUR 


FAC 


6-30 


KN 


(2) 


En  t  ry 

Node  number  to  which  this  load 

1 s  applied : 

CE.l  and  LE.NUKP 


Direction  of  the  load 
component ; 
EQ.2;  Y-translatlon 
EQ.3;  Z-translatlon 

Load  curve  number  that  describes 
the  time  dependence  of  the  load; 
GE.l  and  LE.NLCUR  (Default  -  1) 


Function  multiplier  used  to 
scale  f(t)  for  the  load  at  "i 
EQ.0.0;  default  set  to  "1.0" 


Node  Increment  for  nodal  load 
data  gene  rat  Ion 
EQ.O;  no  generation 


APPENDIX  D.2 
INPUT  DATA  DESCRIPTION 


I 
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|EAD1NC    CARD 

(I)      Begin    each    data    case    with    a    heading    card 

I 

I 

I 

I 

f 

I 

i 

I 

I 

I 

I 

I 

I 

I 

I 
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MASTER  CONTROL  CARDS 
Card  2.1 

1.  The  total  ounber  of  nodes  (NUMNP)  controls  the  amount  of 
data  to  be  read  In  Section  111.   If  NUMNP. EQ.O,  the  program 
terminates  execution. 

2.  The  program  distinguishes  between  linear  and  nonlinear  ele- 
■enta.   Linear  elements  have  their  stiffness  matrices  formed 
only  once,  and  the  formulation  excludes  consideration  of 
either  geometric  or  material  nonl Inea r 1 t les ;  linear  elements 
are  used  to  represent  those  regions  of  a  model  which  can  be 
assumed  to  behave  linearly  and  should  be  used  (where  possi- 
ble) to  Improve  solution  efficiency. 

An  element  group  Is  a  aeries  of  elements  of  a  particular 
type  (e.g.,  TRUSS  or  2-D  CONTINUUM)  in  which  element  numbers 
are  assigned  in  ascending  sequence  beginning  with  "1"  and 
ending  with  the  total  number  of  elements  In  that  particular 
group.   Elements  forming  a  group  must  have  the  same 

a.  Element  type  (NPAR(l) 

b.  Type  of  analysis  (NPAR(3)) 

c.  Material  model  (NPAR(15)) 

d.  Integration  order  (NPAR(IO)) 

Elements  defined  by  the  same  values  of  the  parameters  listed 
above  can  be  broken  down  Into  more  than  one  group.   This 
reduces  solution  efficiency  but  may  be  necessary  to  restrict 
the  size  of  2-D  continuum  element  groups  when  PS-NFAP  is  run 
on  an  IBM-PC. 

It  is  emphasized  that  only  one  type  of  material  model  is 
allowed  for  the  elements  in  a  group  (either  linear  or  non- 
linear).  Also,  it  is  permissible  to  model  a  structure  with 
nonlinear  elements  only,  in  which  case  NEGL.EQ.O  and 
NEGNL.GE.l.   Linear  element  groups  are  Input  before  non- 
linear groups  in  Card  Croup  5. 
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MASTER  CONTROL  CARDS  (Card  2.1  continued) 


13.   Nonlinear  elements  Include  the  effects  of  material  and/or 
geometric  nonl Inearl t les  In  the  formation  of  stiffness  and 
•tress  recovery  Matrices.   The  type  of  nonl lnear 1 t les  to  be 

•    associated  vlth  en  Individual  group  of  elements  Is  defined 
by  Beans  of  data  given  on  the  element  group  control  card. 
In  aeparatlng  elements  Into  nonlinear  groups,  note  that  only 
one  type  of  Material  aodel  Is  allowed  for  the  group.   The 
order  In  which  groups  are  Input  was  discussed  In  Note  (2) 

-    above. 

The  total  number  of  element  groups  processed  by  the  program 
Is  the  sum  NELG+NECNL.   Also,  NEGNL  can  be  zero,  but  then 
the  number  of  linear  element  groups  Bust  be  at  least  one. 

L.   The  MODEX  parameter  determines  whether  the  program  is  to 
check  the  data  without  executing  an  analysis  (lie. 
MODEX. EQ.O)  or  If  the  program  Is  to  perform  an  analysis. 
The  user  has  the  following  choices: 

MODEX  «  0   Data  check;  the  program  will  only  read,  check, 
and  print  input  data. 

1  Perform  analysis  for  a  new  problem. 

2  Perform  analysis  for  a  restart  Job. 

KPER  is  the  number  of  time  periods  which  have  the  same  tl 
step  Incraents,  NPER  Bust  be  at  least  one.   NSTE  Is  the 
number  of  total  solution  (or  load)  steps,  which  Is  calcu- 
lated by  the  program  based  on  TFINAL. 

The  print  Interval  determines  st  which  solution  step  Inter- 
val program  results  are  to  be  printed.   If  IPR1.EQ.4,  output 
Is  produced  at  the  end  of  solution  steps  A,  8,  12,  etc. 

If  IPRI  Is  larger  than  the  total  number  of  solution 
steps  (NSTE),  then  no  output  will  be  printed  during  the 
course  of  solution.   If  IPRI .LE .NSTE ,  then  print  directives 
Bust  be  given  for  displacements  on  Card  2.4  In  this  section. 


me 
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MASTER  CONTROL  CARDS  (Card  2.1  Continued) 

7.  IR1NT  specifies  Che  tine  step  Interval  for  which  the  his- 
tories will  be  ssved  for  restart.   If  1R1NT.EQ.]  the  his- 
tories are  ssved  at  the  end  of  esch  solution  step.   This 
option  is  useful  for  solving  lsrge  scsle  problems. 

8.  ISTOTE  Is  a  parameter  which  specifies  the  maximum  block  size 
to  be  used  for  out-of-core  solution.   If  It  Is  left  blank, 
the  program  will  determine  the  maximum  block  sire.   When 
operating  on  an  IBM-PC  with  640k  of  memory  and  an  A  matrix 
of  size  72010,  ISTOTE  Bust  be  15930  to  limit  the  length  of 
unformatted  records  written  to  disk. 

9.  NRENUM  Indicates  the  node-renumbering  option.  This  minimizes 
the  column  height  of  the  structure  stiffness  thereby  Improv- 
ing solution  efficiency.   The  node  number  16  converted  back 
to  the  original  numbering  scheme  for  the  output.   On  an 
IBM-PC  node  renumbering  1b  done  by  a  separate  program 
(NFMINX)  and  PS-NFAP  reads  the  renumbered  nodes  from  a  file 
named  'NFMIN.NFP'. 
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MASTER  CONTROL  CARDS 
Card  2.2 


1. 


2. 


3  . 


The  stiffness  aatrlx  reforastlon  interval  (ISREF)  Is  ignored 
if  the  aodel  is  composed  of  linear  eleaent  groups  only;  i.e. 
NECNL.EQ.O.   For  linear  problems  the  aatrlx  is  formed  and 
decomposed  only  once. 

For  aodels  containing  nonlinear  element  groups,  the  system 
stiffness  aatrlx  is  reformed  every  ISKEF  solution  iteps 
based  on  conditions  known  at  the  end  of  the  previous  step. 
For  problems  vith  nonlinear  element  groups  using  the  csp 
aaterlal  aodel  ISREF.EQ.l,  i.e.,  the  system  stiffness  aatrlx 
is  reformed  at  the  start  of  every  time  step.   Only  the  non- 
linear portion  of  the  complete  system  stiffness  aatrlx  Is 
reformed;  the  linear  portion  Is  saved  and  reinstated  when 
the  complete  aatrlx  is  calculated. 

Stiffness  reformation  may  also  be  requested  during  equili- 
brium iterations.  Vhen  NUMREF .NE  .0 ,  this  option  is  sctlvated 
for  every  NUMREF  cycles  of  iteration  in  each  time  step. 
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m  dlspla- 
problems 

satlsf ac- 
nstructed 
cap  aodel 


-  D39  - 


MASTER  CONTKOL  CARDS 
Card  2.3 


1. 


2. 


TSPER(I)  is  the  tine  at  which  aolutlon  starts  for  tach  tine 
period.   This  option  provides  the  convenience  to  change  tine 
Step  sire  (DTPER(I))  for  tine  dependent  analysis.  For  a  res- 
tart Job,  the  starting  tine  of  the  first  tine  period  would 
be  the  final  tine  of  the  previous  solution  which  was  alrea 

j   ....  .J   -  _  J  i   -  L 


run  and  saved  on  disk. 


dy 


TFINAL  Is  the  final  time  of  analysis  for  each  run  and  It  Is 
used  to  calculate  the  total  number  of  aolutlon  ateps  (NSTE). 


>STER  CONTROL  CARDS 
rd  2.4 
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1.   For  large  aeihet  It  say  not  be  necessary  to  print  displace- 
•ants  at  every  code.   Hence,  nodes  for  which  printout  la 
dealred  are  grouped  Into  NPB  printout  blocks.   Each  block  of 
Bodes  Is  defined  by  the  node  numbers  of  the  first  and  last 
node  In  the  block  (see  next  card). 

If  NPB.EQ.O  all  nodal  quantities  are  printed  regardless  of 
the  values  of  IDC. 

The  displacement e  at  the  nodes  within  the  blocks  are  printed 
If  IDC.EQ.l. 

If  ghost  elements  are  used  (I.e.,  elements  which  are  Ini- 
tially Inactive),  1DSET  Is  used  to  update  the  Z-direction 
nodal  coordlntes  during  a  large  displacement  analysis. 
IDSET-0  results  In  no  update  of  the  Initial  ghost  element 
geometry.   This  feature  is  useful  for  problems  that  experi- 
ence very  large  settlements. 
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MASTER  CONTROL  CARD 
Card  2.5 


1,   Two  entries  are  expected  for  each  printout  block,  Basely, 
the  first  node  of  the  block  and  the  last  node  of  the  block. 
All  nodal  points  between  these  two  nodes  will  be  Included  In 
the  printout  block.   Maximum  of  6  printout  blocks  are 
allowed . 

Skip  this  card  If  NPB.EQ.O 
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>DAL  POINT  DATA 


rd  3.1 


Modal  data  tutt  be  defined  for  all  (NUMNP)  nodes.   Hode  data 
■ay  be  Input  directly  (I.e.,  each  node  on  lta  own  Individual 


card)  or  the  generation  option 
(aee  note  5).  Admissible  node 
total  number  of  nodes  (NUMNP). 
repeated  or  omitted.  The  last 
"NUMNP". 


■ay  be  used  If  applicable 
numbers  range  from  "1"  to 

Node  numbers  say  not  be 
node  that  Is  Input  Bust  be 


The  print  auppresslon  flag  (PSF)  Is  used  to  eliminate  print- 
ing of  ordered  node  coordinates,  equation  number  assign- 
ments, or  generated  node  coordinates.   The  PSF  character  16 
entered  on  the  first  card  of  nodal  point  data  only. 

Boundary  condition  codes  can  only  be  assigned  the  following 

values  (M  -  2,3) 

ID(M,N)  -  0;  unspecified  (free)  displacement. 

ID(M,N)  ■  1;  deleted  (fixed)  displacement. 

ID(M,N)  -  2;  deleted  (fixed)  displacement,  used  for  generation. 

An  unspecified  [ID(M,N)-0]  degree  of  freedom  16  free  to 
translate  as  the  solution  dictates.   Concentrated  forces  may 
be  applied  In  this  degree  of  freedom. 

One  system  equilibrium  equation  Is  required  for  each 
unspecified  degree  of  freedom  in  the  model.   The  maximum 
number  of  equilibrium  equations  is  always  less  than  two  (2) 
times  the  total  number  of  nodes  in  the  system. 

Constrained  (1D(M,N)-1]  degrees  of  freedom  are  deleted  from 
the  final  set  of  equilibrium  equations.   Deleted  degrees  of 
freedom  are  used  to  define  fixities  (points  of  external 
reaction),  and  any  loads  applied  in  these  degrees  of  freedom 
are  Ignored  by  the  program. 


-  2  can  also  be 
value  on  the  next 


For  generation  purposes  the  value  ID(M,N) 

used.   In  this  case,  If  the  corresponding  value  on  tne  next 
Input  card  Is  aero  (0),  It  is  eet  equal  to  "2".   Considering 
the  deletion  of  degrees  of  freedom  at  two  (2)  have  the  same 
meaning  as  a  plus  one  (+1). 
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NODAL  POINT  DATA 
Card  3.1  (Continued) 


4.   Mode  cards  need  not  be  Input  In  node  order  sequence;  eventu- 
ally, however,  all  nodes  In  the  set  [l.NUMNP]  suit  be 
defined.   Node  data  for  a  aeries  of  nodes 

[N.  .N.  +  l-KNj ,  N1  +  2*KN.,  ....  N£] 

■ay  be  generated  from  information  given  on  two  (2)  cards  In 
eequence-- 


CARD  1  —  Nj ,1D(N1,2),1D(NI,3),Y(N1) 


,KN 


1 


CARD  2  —  N2,ID(N2,2),ID(N2,3),Y(N2),...,KN2 

KN.  Is  the  node  generation  parameter  given  on  the  first  card 
in  the  sequence.  The  first  generated  node  Is  N  +1*KN  ;  the 
second  generated  node  Is  N.+2*KN.  ,  etc.   Generation  contin- 
ues until  node  number  K  -KN   16  established.  Note  that  the 
node  difference  N.-N   Bust  fie  evenly  divisible  by  KN  . 

In  the  generation  the  boundary  condition  codes  (ID(L,J) 
values)  of  the  generated  nodes  are  set  equal  to  those  of 
Node  N  .  The  coordinate  values  (Y,Z)  are  Interpolated 
linearly  . 
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AD  CONTROL  CARDS 
[Card  4.1 


~1 .   KLOAD  determines  the  number  of  cards  to  be  read  as  the  nodal 
forces  In  Card  Group  7.   The  loads  defined  are  conctntratcd 
Bode  forces  that  do  not  change  direction  as  the  structure 
defoms;  i.e.,  the  applied  node  forces  are  conservative 
loads . 
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applied  to  the  structure  by  Beans 
on  [l.e.(  f(t)]  references  and  func- 
d  with  the  loads.   At  tine  t  the 
y  linear  Interpolation  In  the  table 
s  the  aultlpller  is  the  Magnitude  of 
NPTM  is  the  maximum  number  of  (f(t), 
e  any  one  of  the  NLCUR  functions;  an 
have  fewer  than  NPTM  (f(t),  t] 
function  can  be  input  with  more  than 
wo  points  are  required  per  functon; 
In  time  16  not  possible. 


In  elastic-plastic  analysis  with  the  cap  model,  plastic 
loading  and  then  unloading  may  be  encountered.   In  this 
case,  a  switch  from  the  tangent  modulus  to  elastic  material 
matrix  16  necessary  in  order  to  obtain  a  convergent  solu- 
tion. NTURNS  indicates  the  number  of  plastic  unloading  (or 
reverse  loading)  for  the  problem. 
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LOAD  CONTROL  CARDS  (Continued) 
Card  4.2 


1.   TRV  Is  the  list  of  the  time  of  all  the  load  steps  at  which 
the  unloading  begins,  when  the  current  tine  of  the  analysis 
Batch  with  the  tine  indicated  In  TRV  list,  the  elastic 
stiffness  of  the  elements  in  the  specified  group  will  be 
used  in  fording  the  structure  stiffness  latrlx, 
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.USS  ELEMENT  (See  Fig.  V.l) 
>rd  5.1 


1.   TRUSS  element  numbers  begin  with  one  (1)  and  end  with  the 
total  number  of  elements  In  this  group,  NPAR(2). 

The  parameter  KPAR(3)  Is  applicable  only  If  the  element 
group  Is  nonlinear.  If  NPAR( 3 ) . EQ. 1 ,  no  geometric  nonlinear- 
lties  are  taken  Into  account,  I.e.  the  geometric  atlffness 
matrix  is  not  Included.   If  NPAR(3).EQ.2  large  displacement 
effects  are  Included  In  the  analysis,  but  small  strains  are 
assumed  In  the  calculation  of  clement  forces. 

In  any  one  element  group  only  one  material  model  can  be 
used,  and  this  model  type  Is  defined  by  the  entry  NPAR(15). 
If  NPAR(15).EQ.l  the  model  Is  defined  by  Young's  modulus 
only  and  NPAR( 1 7 ) . EQ.O.   If  NPAR( 1 5) .EQ.2  the  stress-strain 
curve  is  defined  by  Input  data. 

The  model  defined  for  the  element  group  must  be  con- 
sistent with  the  nonlinear  formulation  used  (defined  by 
KPAR(3))  and  the  requirement  of  equilibrium  iteration  as 
defined  on  Card  2.2  of  the  Master  Contrl  Cards  (Card  Croup 
2).  As  stated  In  note  (3)  of  card  2.2,  equilibrium  Itera- 
tions can  only  be  performed  if  the  model  allows  for  itera- 
tion, and  if  at  least  one  nonlinear  element  group  is  used  in 
the  analysis . 

The  variable  NPAR(16)  defines  the  number  of  sets  of 
material/section  properties  to  be  read  in. 

For  the  nonlinear  elastic  material  model  (NPAR( 1 5) .EQ. 2 ) , 
MPAR(17)  Is  two  times  the  maximum  number  of  points  used  to 
describe  the  nonlinear  stress  strain  curves. 
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♦  Sign:  TentHe  Force 


Figure   V. 1 .      A  truss   clement. 
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RL'SS  ELEMENT  (Continued) 
ard  5.2 


|l .   MPAR(lfc)  different  linear  elastic  saterlals  are  input  In 
this  section,  provided  NPAR( 1 5 ) . EQ. 1 .   Note  that  one 

Iaaterlal/sectlon  la  defined  to  have  the  line  Young's 
aodulus,  area,  and  Initial  axial  force. 


I 
f 
I 

I 
f 
I 

i 

i 
f 
I 
I 
I 
I 
I 
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TRUSS  ELEMENT  (Continued) 
Card  5.3 


1. 

2. 
3. 


One  section  property  card  Is  defined  to  have  the  lice  area, 
Intlal  axial  force,  and  tine  when  element  Is  activated. 

P1N1T  Is  an  Initial  axial  force. 

Tine  element  Is  activated.   This  option  Is  useful  for  Incre- 
mental construction. 

The  stress-strain  curve  is  defined  by  straight  lines  between 
the  Input  points  (c  ,  o  ).   From  the  at ress-st rain  curve 
total  stresses  and  the  tangent  modulus  are  evaluated  for  a 
given  atraln  (see  Fig.  V.2). 

The  variable  NCON  was  defined  in  CARD  5.1  by  the  variable 
NPAR(17). 

This  model  can  only  be  used  In  a  nonlinear  element  group. 


I 
I 
1 
1 
I 
I 
I 
I 
I 
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Figure  V.2.   A  nonlinear  »tress-»traln  curve  for  tne  trus« 
e 1 ement . 
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TRUSS  ELEMENT 
Card  5.4 


1.   tcfer  to  Figure  V.l 
•ode  Is  eero. 


For  two-node  trues  eeitberi  the  third 


2*   Elements  Bust  be  Input  In  Increasing  elesent  number  order. 
If  cards  for  elesents  (M+l,  M* 2 , . . .  ,M« J ]  are  omitted,  these 
"J"  alsslng  elenenta  are  generated  using  MTYP  of  elesent  "M" 
•nd  by  Incrementing  the  node  numbers  of  successive  elements 
with  the  value  "KG";  KG  la  taken  from  the  first  card  of  the 
element  generation  aequence  (I.e.,  from  the  "M-th"  element 
card) . 


2-E>  CONTINUUM  CLEMENTS 

|-D  CONTINUUM  elenentt  arc  4-  to  B-node  ieoparanet ri c  plane 
•train  quadrilaterals.   They  amet  be  input  in  the  global  Y-2 
^lane.   Bach  eleacnt  node  hat  two  (2)  tranalat lonal  degree*  of 
flreedoa. 

• 
I 

1 
I 
I 

i 


P 
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2-D  CONTINUUM  ELEMENTS 
Card  6.1 


(1)  2-D  CONTINUUM  element  nuibert  begin  with  mat  (1)  and 
•nd  with  total  number  of  elements  In  this  group, 
*PAR(2).   Element  data  arc  Input  in  CARD  »•*•   For  the 
11M-PC  version  NPAR(3)*)56  ♦  NPAR( 1 6)*46 .LE.l S»36  to 
Halt  the  length  of  unforaatted  records   written  to 
disk.   Larger  groups  aay  have  to  be  divided  into  2 
groups. 

(2)  NPAR(3)  Is  applicable  for  non-linear  groups  only  and 
determines  If  geoaetrlcal  nonllnearl t las  are  to  be 
Included  In  the  analysis.   If  NPAR( 3) .EQ.l  displa- 
cements and  strains  are  assumed  to  be  Infinitesimal. 
In  the  updsted  Lagranglan  Formulation  all  geometric 
effects  (large  displacements  or  large  strains)  are 
Included  In  the  analysis. 

(3)  NPAR(7)  limits  the  number  of  nodes  that  can  be  used  to 
describe  any  of  the  elements  In  this  group.  A  minimus 
of  4  and  a  maximum  of  6  nodes  are  used  to  describe  the 
2-D  CONTINUUM  elements.  Although  5  or  7  node  elements 
can  be  used,  to  ensure  numerical  accuracy  either  4,  6 
or  6  node  elements  are  recommended  for  use  as  ahovn  In 
Fig.  VI. 1  . 

(6)   For  rectangular  elements  an  Integration  order  of  "2"  Is 
sufficient.   If  the  element  Is  distorted,  a  higher 
Integration  order  need  be  used.  Notice  that  apart  from 
the  larger  computational  effort  in  the  calculation  of 
the  element  matrices,  more  working  storage  is  required 
if  a  nonlinear  material  model  Is  used  (see  note  (5) 
below).   Stresses  are  output  at  each  Integration  loca- 
tion as  shown  in  Fig.  VI. 2. 

(5)   Only  one  asterlal  model  (defined  by  the  value  of 

NPAR()5))  Is  allowed  in  an  element  group.  PS-NFAP  sup- 
ports only  model  numbers  1  and  10.   If  NPAR(15)  Is  Ml" 
the  entries  for  NPAR(17)  and  NPAR(IB)  are  Ignored  by 
the  program.   If  NPAR(15)  Is  "10",    material  constants 
are  read  Into  the  storage  array  "PROP"  which  has  dimen- 
sions NPARU7)  by  NPAR(16):  I.e.,  property  constants 
are  atored  as, 

((PR0P(1,J)(  X  -  1,  IIP Al(  17)),  J  -  1,  NPAR(U)) 

Model  10  requires  that  NPAR( 1 7) .EQ.22  and  that  the 
dimension  of  the  storage  array  for  element  history 
NPAR( )6).EQ.12.  The  NPAR(IB)  parameters  which  charac- 
terise element  hiatory  suit  be  retained  for  materially 


-  VbU    - 


2-D  CONTINUUM  ELEMENTS 
Card  6.1  (Continued) 


nonlinear  eleaenta  to  that  properties  can  be  chosen  In 
the  current  solution  step.  The  working  storage  array 
(WA)  is  dimensioned  by  1DWA*NPAR( 2 ) ;  where: 

1DWA  -  NPAR(18)*(NPAR(10)**2)*NPAR(2)*2 


The  model  defined  for  the  element  group  oust  be  con- 
sistent with  the  nonlinear  formulation  used  (defined  by 
NPAR(3)),  the  element  type  (defined  by  NPAR(5)),  and 
the  requirement  of  equilibrium  Iteration  as  defined  on 
Card  2.2  of  the  Master  Control  Cards.  As  stated  In  note 
(3)  of  Card  2.2,  equilibrium  iteratons  can  only  be  per- 
formed if  the  model  allows  for  iterations,  and  if  at 
least  one  nonlinear  element  group  is  used  in  the 
analysis . 
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4-nod#  •lament 


6      6-nod«  sltmtnt 


6  6-nod*  slsmtnt 


Figure  VI.  I.      Recommended  use  of  2-D  elements 
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Strtss  Printout  Tsbls  Points 
for  2-D  Continuum  Elsmsnts 


2-Polnt  Inttgntlon 


3-Polnt  Inttgntlon 


Figure  VI. 2.   Stress  print-out  convention  for  2-D  continuum 
elements. 
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2-D  CONTINUUM  ELEMENT 


Card  6.3 
(1) 


1 


MODEL  1  Is  a  linear  clastic,  Isotropic  Material  defined 
by  two  (2)  positive  constants  (E,v);  I.e.,  If 

NPAR( J  5). Eg.] ,  NPAR(17)  Is  let  to  "2"  by  default. 
MODEL  1  can  be  used  with  linear  or  nonlinear  eleeent 
groups.  Since  the  material  constants  are  Independent  of 
history,  NPAR(18)  Is  set  to  "0"  by  default. 
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•D  CONTINUUM  ELEMENT 

:ard  6.3.10   Soil  Plasticity  Model,  NPAK(15)  -  10 

node]  10  Is  a  capped  plasticity  model  consisting  of  a  work- 
hardening  aurface  and  a  perfectly  plastic  surface.   The  work- 
Hardening  surface  has  the  form: 

F  -  (lj  -  XL)2  +  R2J2  -  (X  -  XL)2  -  0 

jfhere  I.  la  the  first  stress  lnvarlent  and  J.  is  the  second  devi- 
atorlc  stress  lnvarlent.  The  perfectly  plastic  surface  is  of  the 
^rucktr-Pnger  type  and  has  the  form: 

*  AM  *  lj  4  J*/2  -  AC  -  0 

■'he  cap  model  is  shown  in  Fig.  VI. 3. 

(1) 

*The  parameters  K, ,  K„,  A  ,  K  .  ,  C.  and  G.  are  used  to 
1    2.         pmln    1       2 


represent  the  elastic  moduli.   They  have  the  form: 

•A  *[-!  /(3*A  )1   ,  . 
1   p     1      p        Bin 


K  -  MAX  (K  *A  *[-I./(3*A  )]   ,  K_,_) 


G  -  C   +  G   *  K 

where  K  and  G  are  the  bulk  and  shear  moduli,  respec- 
tively and  K.  ,  K_,  G.,  and  G.  are  fitting  parameters; 
K  .   must  be  .GE.0.0  and  is  the  minimum  value  of  K;  A 
is  atmospheric  pressure  and  has  the  same  units  as  K. 

(2)  AM  and  AC  are  defined  In  Figure  VI. 3. 

(3)  AW  and  AD  are  material  parameters  for  the  hardening 
rule. 

Ckk  "  AW  (exP(AD*x>  "  ?•] 

(4)  R  is  the  aspect  ratio  of  the  work-hardening  elliptical 
cap.  The  aspect  ratio  Is  defined  as  the  ratio  of  the 
horizontal  to  vertical  axis  of  the  ellipse. 


(5)  XL  Is  the  initial  position  of  the  elliptical  cap 
defined  by  the  intersection  of  the  cap  and  ultimate 
failures  aurface.   For  normally  consolidated  soil  XL 
should  be  aet  equal  to  41.  and  the  initial  cap  position 
will  be  calculated  by  the  program. 

(6)  TENCUT  la  the  maximum  principle  tension  stress.   The 
material  la  aasumed  to  fracture  after  reaching  thla 
value. 
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2-D  CONTINUUM  ELEMENT 
Card  6.3.10  (continued) 

(7)  STATE  Is  a  Material  parameter  which  controls  whether 
the  cap  Is  allowed  to  contract  as  well  as  expand  along 
the  hydrostatic  axis: 

STATE  «  1  for  expansion  and  contraction. 
>  1  for  expansion  only 

(8)  A.,  A   and  K   are  used  to  define  the  Initial  stresses: 

o   -  A*Z  4  A, 

z  1       2 

o   -  K  *o 
y     0   * 

o   -  K  *o 
x     0   * 

where  Z  Is  the  Z-coordlnate  of  the  Gauss  point  under 
consideration.   If  Initial  stresses  are  specified,  they 
oust  be  compensated  by  balancing  nodal  loads  In  Section 
VII.   Initial  stresses  cannot  be  specified  If  Incremen- 
tal construction  (PR0P(20),  PR0P(21),  and  PROP(22))  Is 
used . 

(9)  SWITCH  is  the  time  at  which  the  analysis  Is  switched 
from  drained  to  undrained  conditions.   It  must 
correspond  to  the  exact  time  at  which  one  of  the  solu- 
tion time  steps  occurs.   For  times  less  than  switch  no 
pore  pressures  are  generated.   Set  SWITCH .EQ.O .0  If 
this  option  Is  not  used. 


(10)  FAC  Is  used  to  compute  the  bulk 
use  In  an  undrained  analysis: 


lodulus  of  water  for 


K   -  FAC*K 

w 

where  K  is  the  bulk  modulus  of  the  soil.   The  pore 
pressures  are  assumed  to  be  0  for  material  that  has 
failed  In  tension  so  the  program  sets  FAC  -  0  for  this 
case. 


FAC  >  0   computes  positive  and  negative  excess  pore 

p  ressures 
FAC  -  0   drained  conditions;  no  pore  pressures  generated 
FAC  <  0   computes  only  negative  excess  port  pressures; 

sets  pore  pressure  -  0  for  material  vlth  positive 

volumetric  strain 
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2-D  CONTINUUM  ELEMENT 
Card  6.3.10  (continued) 


(11) 


(12) 


TIME-ON  Is 
vlth  th: 
TIME-ON 


the  tine  at  which  the  elements  associated 
vlth  this  set  of  aaterlal  parameters  becomes  active. 

just  correspond  to  the  exact  time  at  which  one 
of  the  program  time  steps  occurs.   For  time  less  than 
TIME-ON,  no  element  stresses  develop.   This  is  used  for 
Incremental  construction.   Incremental  construction 
cannot  be  used  If  initial  stresses  (PR0P(15),  PR0P(16), 
and  PR0P(17))  are  specified. 

WGT  and  BOUY  are  used  to  calculate  the  gravity  and 
bouyancy  forces  on  an  element  for  which  the  TIME-ON 
option  Is  used.   Bouyancy  forces  are  calculated  if  a 
Gauss  point  falls  below  the  water  table  which  Is 
assumed  to  be  at  Z  ■  0.0.   The  weight  and  bouyancy 
forces  are  Increased  from  zero  at  time  *T1ME-0N'  to  the 
value  defined  by  (WGT-BOUY)  at  time  "TIME-ON"  +  1  in 
the  manner 

*   BF  -  XFAC  *  (WGT-BOUY) 

XFAC  -  Time-(TIME-ON) 

XFAC  -  Minimum  of  (XFAC,  1.0) 


Only  two  cases  can  be  modeled:  (1)  the  water  table  at  the  ground 
aurface  In  which  case  BOUY  is  the  unit  weight  of  water;  or  (2)  no 
water  table  in  which  case  BOUY  Is  zero. 

Note:   For  this  Model   NPAR( 1 7) .EQ.2 2  and  NPAR( 1 8) .EQ. 1 2 
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!-D    CONTINUUM    ELEMENT    (Continued) 
lard    6.4 


(1) 


(2) 


Elements  luit  be  Input  In  ascending  element  ouaber 
order.   If  data  cards  for  elements  [M+l,  M+2,  •••»  M+J] 
are  omitted,  these  "J"  alsslng  elements  are  generated 
using  1EL,  IPS,  and  MTYP  given  on  the  card  for  element 
"M"  and  by  incrementing  node  numbers  of  successive  ele- 
tenti  vith  the  value  "KG";  the  value  of  KG  used  for 
i cerement  at  ion  Is  taken  from  M-th  element  card,  and 
only  the  non-zero  nodes  appearing  on  the  M-th  element 
card  are  lncreaented  when  generating  Biasing  element 


data. 


An    example    of    a    6-node    element    (1EL.EQ.6)    shown    below: 


Noda  Location  5 

Not  Used 


Nodt  Location  7 
Not  Ustd 
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Nod* 

r 

• 

1 

♦1 

♦  1 

2 

-1 

♦  1 

3 

-1 

-1 

4 

♦1 

-1 

5 

0 

♦1 

6 

-1 

0 

7 

0 

-1 

8 

♦  1 

0 

Figure  VI .4 


Element  node  number  Input  sequence  for  2-D 
continuum  elements. 
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LOCATIONS  1,  2.  5 


Figure  VI .5. 


Element  node  numbering  for  triangular 
elements. 
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APPLIED  LOADS  DATA 
Card  7.1b 


(1)  Time  values  at  luccettlve  point*  auat 
Increase  In  aagnltude  (I.e.,  T1MV(1)  < 
T1MV(2)  <  TIMV(3),  etc.).   The  last  else 
value  Bust  be  greater  than  or  equal  to  the 
tine  at  the  end  of  solution;  I.e.  TIMV(NPTS) 
>  TSTART  ♦  NSTE*DT 

(2)  Input  as  aany  cards  In  this  section  as  are 
required  to  define  NPTS  points,  four  points 
per  card. 
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APPLIED  LOADS  DATA 
Card  7.2 


(1)  If  the  same  degree  of  freedom  (IDIRN)  at  the 

eaae  code  (NOD)  la  given  a  multiple  number  of 

tines,  the  program  adds  the  loads  algebrai- 
cally. 

(2)  Nodal  loads  data  which  are  applied  as  a 
aeries  of  codes 

»j,  Nj  4  1*  KN,  Nj  ♦  2*KN ,  N 

■ay  be  generated  from  Information  given  as  2 
cards  In  aequence. 

Card  1   /H.,  IDIRN,  NCUR,  FAC,  KN  / 

Card  2   /N2 ,  IDIRN,  NCUR,  FAC/ 

KN  Is  the  generation  parameter  given  on  the 
first  card.   The  first  generated  loaded  node 
Is  N   +  1*KN;  the  second  generated  loaded 
node  Is  N   *  2*KN,  etc.   Generation  continues 
until  node  number  K  -KN  Is  established.   The 
node  difference  N  -K   must  be  evenly  divisi- 
ble by  KN.        l       l 


The  aeries  of  generaed  nodes  will  have 
the  same  IDRIN,  NCUR,  and  FAC. 


APPENDIX  E 

EXAMPLE 

ANALYSIS  OF  REINFORCED  EMBANKMENT 

ON  SOFT  GROUND 
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INTRODUCTION 

This  appendix  Illustrates  the  use  of  the  computer  pro- 
grams described  In  this  manual  and  procedures  used  to  ana- 
lyze reinforced  embankments  constructed  on  soft  foundations. 
Example  Input  and  output  files  for  programs  CPCALC*  CAP, 
AUTOGEN,  and  PS-NFAP  along  with  Its  auxiliary  program  NFMINX 
•re  given.   Appendices  A  through  0  should  be  consulted  for 
detailed  explanation  of  the  Individual  records  shown  In  the 
example  Input  files.   A  detailed  description  of  the  proce- 
dures to  determine  the  cap  parameters  Is  given  In  Humphrey 
and  Holtz  (1986)i  only  parameters  that  are  directly  related 
to  program  operation  are  discussed  here. 

PROBLEM  DEFINITION 

The  problem  chosen  for  this  example  Is  a  granular 
embankment  with  a  90-foot  base  width  and  2htlv  side  slopes. 
The  foundation  Is  composed  of  a  7. 5-foot  thick  overconsol I- 
dated  dried  crust  underlain  by  22.5  feet  of  soft  normally 
consolidated  clay  as  shown  In  Fig.  E.l.   The  water  table  1s 
at  the  ground  surface.   The  embankment  was  constructed  In 
horizontal  lifts  and  there  was  no  time  for  excess  pore  water 
pressures  In  the  foundation  to  dissipate.   The  fill  has  ae 
unit  weight  of  125  psf.  Mohr-Coulomb  friction  angle  of  32  • 
and  negligible  cohesion.   Elastic  parameters  were  based  on 
nonlinear  stress  strain  parameters  given  by  Rowe,  et  el. 
(1984).   A  Polsson's  ratio  of  0.33  was  used.   Properties  for 
the  normally  consolidated  foundation  soils  ere  given  1n 
Table  E.l.   The  strength  parameters  were  obtained  from  K 
consolidated  undralned  trlexlel  compression  tests.   The  the 
surface  crust  has  the  same  properties  except  that  the 
undralned  shear  strength  (s  )  Is  125  psf  and  K   Is  0.6.   The 
embankment  Is  reinforced  byua  single  layer  of  tensile  rein- 
forcement placed  on  the  ground  surface  as  shown  In  Fig.  E.l. 
The  reinforcement  has  a  modulus  of  60  k/ft  (5000  lb/In.) 
which  Is  representative  of  a  strong  geotextlle  or  geogrld. 

CAP  PARAMETERS  FOR  FOUNDATION 

Behavior  of  the  foundation  soil  was  represented  with 
the  cap  model.   To  determine  o  from  ♦'  and  J2f   'cvo 
from  s   It  Is  necessary  to  assume  a  value  of  r\f   which  is  a 
functl8n  of  the  Intermediate  principal  stress  at  failure 
(Humphrey  end  Holtz,  1986)  where  rjf  I*  defined  as 

nf  -  (o2  -  oi)f  I   (c3  -  o,)f 

wheres     o   ■  major  principal  stress 

o«  »  intermediate  principal  stress 

Oj  ■  minor  principal  stress 
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The  assumed  nf  Is  compared  to  the  average  value  at  failure 

In  the  foundation  as  computed  by  PS-NFAP.   The  resulting  o  * 

0.183  end  J  ''  fa'      «  0.339. 
if  vo 

The  cap  parameters  were  determined  with  the  aid  of  the 
program  CPCALC.  The  program  calculated  D*  ,  W/(a-b),  and  R 
as  well  as  x./x    ,  R__  .  and  (x-/x  )    .   pPogram  Input  and 

output  Is  shown  in  Fig.  E.2  end  the  ver lebles  ere  defined  In 
Append \x   A . 

The  remaining  cap  parameters  were  determined  using  the 
procedures  given  1n  Humphrey  and  Holtz  (1966)  and  the  guide- 
lines given  1n  Appendix  D.   The  parameters  for  the  normally 
consolidated  end  overconsol Idated  soils  ere  summarized  In 
Table  E.2.   The  velue  of  D  Is  e  function  of  x      which  can  be 
computed  once  the  FE  mesh  Is  selected.   Note  €hat  a  minimum 
value  of  the  bulk  modulus  (K    )  of  1.0  ksf  was  assumed. 
This  small  value  avoids  numerical  difficulties  et  low  stress 
levels.   A  small  cohesion  («)  Is  el  so  used  to  avoid  numeri- 
cal difficulties.   The  water  table  is  at  the  ground  surface 
so  v  Is  the  buoyant  unit  weight.    The  Initial  stresses  ere 
specified  so  balancing  nodal  loads  must  be  applied  as  dis- 
cussed In  a  leter  section.   The  options  to  switch  from 
dretned  to  undrained  analysis  and  for  Incremental  construc- 
tion ere  not  used  so  SWITCH,  T1ME0N,  WGT,  end  BOUY  *  0.0. 

The  undreined  stress  strain  behavior  predicted  using 
the  cap  parameters  In  Table  E.2  was  checked  with  the  program 
CAP  for  a  normally  consolidated  sample  with  an  Initial  ver- 
tical stress  equivalent  to  a  depth  of  10  feet  (A_  ■   o'   « 
0.526  ksf)  and  e  K   of  0.47.   The  Input  file  Is  Shown  YR 
Fig.  E.3.   Appendix  B  should  be  consulted  for  a  description 
of  the  input  records.   The  output  file  Is  listed  In  Appendix 
E.l  elong  with  2  files  containing  stress  strain  data  for 
Input  Into  commercial  plotting  programs.   The  calculated 


J»£  /P. 


and  normalized  pore  pressure  versus  axial 
1s  plotted  In  Fig.  E.4.   The  calculated 


w-f  'wvo  was  0-340  which  1s  within  0.3%  of  the  value 
used  to  calibrate  the  model.   These  should  be  compared  to 
observed  laboratory  test  results  and  the  model  parameters 
adjusted  If  necessary  to  obtain  a  better  fit. 

PARAMETERS  FOR  EMBANKMENT  FILL 

The  embankment  fill  was  modeled  using  only  the  Drucker- 
Preger  ultimate  failure  surface  (Humphrey  and  Holtz,  1986). 
This  Is  done  by  fixing  the  Initial  position  of  the  cep  (XL) 
well  beyond  the  maximum  stress  level  that  occurs  In  the 
problem.   An  XL  of  -1.0x10    ksf  was  chosen  and  STATE  ■  2.0 
prevents  contraction  of  the  cap.   The  remaining  parameters 
directly  associated  with  the  cap  (AW,  AD,  end  R)  are  not 
used  and  are  assigned  e  value  of  1.0.   Other  parameters  were 
determined  from  the  soil  properties  given  above  using  the 
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procedures  given  In  Humphrey  and  Holtz  (1986)  and  from  the 
guidelines  given  fn  Appendix  D.   All  the  parameters  are  sum- 
marized In  Table  E.3.   Small  values  of  K  .   and  c  are  as- 
sumed to  avoid  numerical  difficulties,   incremental  con- 
struction was  used  so  WGT  and  BOUY  mre   specified.   Selection 
Of  TIHEON  Is  discussed  In  the  following  section. 

SELECTION  Of  FINITE  ELEHENT  HESH 

The  problem  geometry  was  modeled  by  the  finite  element 
Mesh  shown  In  Fig.  E.5.   The  geometry  Is  symmetric  about  the 
embankment  center  line  so  only  half  the  problem  needs  to  be 
represented  by  the  FE  mesh.   The  ground  surface  should  be 
chosen  as  Z  »  0.0.   Two  nonlinear  element  groups  were  used 
so  (NEGNL  «  2  on  PS-NFAP  card  2. It  see  Appendix  D.l).   One 
modeled  the  foundation  and  embankment  fill  using  52  nonlin- 
ear 2-D  elements  and  the  other  modeled  the  reinforcement 
using  3  nonlinear  TRUSS  elements.   There  mre   no  linear 
element  groups  (NEGL  »  0  on  PS-NFAP  card  2.1). 

The  foundation  was  divided  into  4  layers  of  6  node  non- 
linear 2-D  elements  with  dimensions  7.5  feet  high  by  15  feet 
long.   The  top  layer  represents  the  dried  surface  crust  and 
the  bottom  3  layers  represents  the  normally  consolidated 
soils.   The  value  of  x   can  now  be  computed  using  the  proce- 
dures given  in  Humphre?  and  Holtz  (1986)  and  then  used  to 
determine  the  cap  parameter  D  In  each  layer. 

The  embankment  was  constructed  In  horizontal  lifts. 
This  was  simulated  in  the  FE  analysts  with  an  Incremental 
procedure.   Each  lift  Is  represented  by  a  row  of  elements. 
During  each  time  step  In  the  analysts  the  stresses  In  the 
next  row  of  elements  mre    Increased  from  zero  to  their  full 
value  In  one  or  more  sub- Increments.   Elements  above  the 
current  lift  mre    Inactive  and  mre   said  to  be  "ghost  ele- 
ment s".   The  cap  parameter  TIHEON  specifies  the  time  when  an 
element  Is  activated.   For  example  the  stresses  In  the  first 
lift  are  increased  from  zero  at  time  »  0.0  to  their  full 
value  at  time  »  1.0  so  the  elements  have  TIHEON  •  0.0  s  o.oi 
(note  that  a  value  slightly  greater  than  zero  should  be 
specified  for  the  first  lift  to  prevent  It  from  being 
Included  In  the  Initial  conditions  which  are  established  at 
time  •  0.0),  for  the  second  lift  TIHEON  ■  1.0  and  the 
•tresses  mre   fully  applied  at  time  -  2.0,  for  the  third  lift 
TIHEON  ■  2.0,  etc.   A  lift  thickness  of  1.B75  feet  was  used 
and  the  load  was  Increased  to  its  full  value  in  2  sub- 
Increments  (DTPER  -  0.5  on  PS-NFAP  card  2.3).   Although  the 
thickness  Is  larger  than  typical  construction  lifts*  thinner 
lifts  would  require  either  unacceptably  thin  elements  or  an 
excessive  number  of  elements.   It  was  estimated  that  the 
embankment  would  fall  at  a  height  of  less  than  11.25  feet  so 
6  lifts  were  used.   A  combination  of  8  node  rectangular, 
trapezoidal  and  triangular  nonlinear  2-D  elements  were  used 
as  shown  on  Fig.  E.5. 
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The  reinforcement  was  modeled  with  three  three-node 
nonlinear  TRUSS  elements  located  at  the  Interface  between 
the  fill  and  foundation  tolls  as  shown  on  Fig-  E. 5.   To  sim- 
ulate a  material  that  can  support  only  tensile  forces  the 
•tress  strain  curve  shown  In  Fig.  E.6  was  used.   Stress  Is 
expressed  as  the  force  times  cross- sect  I onal  mrmm.      To  be 
consistent  the  cross-sectional  area  of  the  elements  fs  1.0 
sq.ft.  (AREA  ■  I.Don  PS-NFAP  card  5.3a).   Note  that  the 
Curve  does  not  pass  exactly  through  the  origin  since  this 
can  result  In  numerical  difficulties  under  some  circum- 
stances.  The  reinforcement  Is  In  place  at  the  Start  of  con- 
struction SO  TIHEON  »  0.0  (TIHEON  «  0.0  Oh  PS-NFAP  card 
5.3a). 

USE  AUTOGEN  TO  CREATE  PS-NFAP  INPUT  FILE 

AUTOGEN  was  used  to  create  the  PS-NFAP  Input  file.   The 
Input  and  output  flies  mre   shown  In  Appendices  E.2  and  E.3 
respectively.   Comments  explaining  how  the  AUTOGEN  commands 
were  used  mre    Included  on  the  Input  file.   Appendix  C  should 
be  consulted  for  a  detailed  description  of  the  commands  and 
PS-NFAP  input  records  mre   described  In  Appendix  D.   The 
output  file  will  be  used  as  the  PS-NFAP  Input  file. 

The  LOAD  command  was  used  to  generate  nodal  loads  to 
balance  the  Initial  vertical  stress  specified  In  the  founda- 
tion soil  by  the  cap  parameters.   A  vertical  load  1s  gener- 
ated for  each  of  the  136  nodes  In  the  foundation  (therefore 
NLOAD  «  136  on  PS-NFAP  card  4.1).   These  loads  correspond  to 
load  curve  1  (PS-NFAP  cards  7.1a  and  7.1b)  as  shown  In  the 
AUTOGEN  Input  file. 

Load  curves  2   through  7  mre   optional  and  serve  only  as 
a  reminder  of  the  load  applied  In  each  construction  lift. 
For  example*  curve  2   corresponds  to  the  first  lift  and  shows 
that  the  load  Is  zero  at  time  ■  0.0  and  Is  -1.0  (the  minus 
sign  indicates  that  the  force  is  applied  downward)  for  time 
2.  1.0.   The  last  card  group  In  the  AUTOGEN  Input  file  Is  the 
FOLL  command. 

USE  NFMINX  TO  RENUMBER  NODES 

The  auxiliary  program  NFMINX  was  used  to  renumber  the 
nodes  to  minimize  the  bandwidth  of  the  stiffness  matrix. 
This  separate  step  Is  only  necessary  In  the  IBM-PC  version; 
It  Is  done  directly  by  PS-NFAP  In  the  mainframe  version. 
The  Input  file  Is  the  PS-NFAP  Input  file  shown  In  Appendix 
E.3.   HFMINX  generates  a  list  of  the  renumbered  nodes  shown 
In  Appendix  E.4  and  a  file  named  *NFMIN.NFP'  which  will  be 
used  by  PS-NFAP. 
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EXAMPLE  RUN  WITH  PS-NFAP 

The  problem  Is  now  ready  for  analysts  with  PS-NFAP. 
The  so  Jut  Jon  was  carried  out  up  to  time  ■  1.0  during  which 
the  first  lift  was  placed.   The  Input  file  Is  shown  In 
Appendix  E.3  and  the  output  file  Is  shown  In  Appendix  E.5. 
At  the  end  of  program  execution  the  files  TAPES  and  TAPE9 
•re  saved  on  disk  for  future  use  by  a  restart  Job  that  would 
resume  the  analysts  at  time  «  1.0.   If  the  analysts  had  been 
continued  It  would  show  that  the  solution  does  not  converge 
at  time  ■  5.5  so  5  lifts  can  be  placed  prior  to  embankment 
fa  1 1  ur  e . 

The  problem  was  run  on  an  IBM-PC  with  with  a  math  co- 
processor chip,  640K  of  memory,  and  a  hard  disk.   The  sire 
of  the  A-matrlx  In  PS-NFAP  was  set  at  72010.   The  program 
used  the  out  of  core  solution  option  and  the  execution  time 
was  about  20  minutes. 

Comments  were  added  to  the  output  shown  In  Appendix  E.5 
to  assist  In  Its  Interpretation.  In  addition,  the  following 
should  be  kept  in  mind: 

1.  The  coordinate  system  Is  Y-dt recti on  horizontal 
(positive  to  right),  2-dlrectlon  vertical  (negative  Is 
down)  ,  and  X-dt recti  on  parallel  to  the  longitudinal 
axts. 

2.  Compressive  stresses  and  strains  are  negative.   This  Is 
opposite  to  the  convention  used  In  soil  mechanics. 

3.  Stresses  are  output  as  effective  stresses. 

Equilibrium  Iteration  Output 

PS-NFAP  outputs  equilibrium  iteration  Information  for 
each  time  step  (see  Appendix  E.5).   The  variables  are 
defined  as  follows: 

LOAD  VECTOR  NORM  «  [Z (external ly  applied  nodal  loads)2)''2 

INCREMENTAL  LOAD  NORM-  ri( loads  added  during  time  step)2) 

MAX  INCREMENTAL  LOAD  -  Maximum  of  all  Incremental  loads 

added  at  the  nodes 

SUM  GENERATED  GRAV  LOAD  «  Sum  of  total  load  applied  using 

Incremental  construction 

UNBALANCED  NORM  ■  measure  of  unbalanced  forces        9    . /? 

»  [Kapplled  load  -  resisting  forces)  3 
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MAX  UNBALANCED  —  DOF  -  maximum  unbalanced  force 

•  MAX(appl led  load  -  resisting  force) 
and  the  DOF  or  equation  number  at 
which  It  occurs 

IDNORM  «  Incremental  displacement  norm,  j  .^ 
■  IK  Incremental  displacements)  ) 

DNORM  «  Displacement  norm       -  . /2 
-  (Knodal  displacements)  J 

GTOT  ■  Sum  of  total  load  applied  using  Incremental 
construction  minus  buoyancy  forces 

If  the  solution  Is  converging  properly  the  UNBALANCED  NORM, 
magnitude  of  the  MAX  UNBALANCED  force,  and  IDNORM  should 
decrease  In  each  successive  Iteration  and  DNORM  should 
approach  a  constant  value. 

Stress  Output  for  2-D  Elements 

Stresses  are   output  at  the  iteration  points  as  shown  in 
Appendix  D.  Fig.  VI. 2.   The  variables  In  the  STRESS 
CALCULATIONS  FOR  ELEMENT  GROUP  1  (2/D  CONTINUUM)  are  as 
fol lows: 

NUM        ■  Element  number 

IPT        «  Iteration  point  number 

STRESS  STATE  ■  State  of  stress;  on  cap,  elastic,  etc. 

STRESS-XX  «  o'   •  Effective  normal  stress  In  X-dlrectlon 

XX 

STRESS-YY   -  o'   «  Effective  normal  stress  In  Y-dl recti  on 
yy 

STRESS-Z2   ■  o'   »  Effective  normal  stress  In  2-dlrectlon 

STRESS-Y2   «  o'   ■  o'   ■  Shear  stress  on  Y  plane  In 
**     zy    2-dlrectlon 

PORE  P     «  Excess  pore  water  pressure 

MAX  STRESS  «  o'    «  Major  effective  principal  stress 

MIN  STRESS  -  o'  -  Minor  effective  principal  stress 

ANGLE      «  Orientation  of  principal  axes  relative  to  Y  and 
2  axes  In  degrees;  counterclockwise  Is  positive 

ETA        ■  n  •  (o2   -  Oj)  /  <o3  -  0j) 

CAP  POSITION  -  x  »  Intersection  of  cap  with  I J  axis 
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Stress  output  for  TRUSS  elements 

The  force  In  the  reinforcement  Is  output  at  3  Integra- 
tion points  (IP)  In  each  element.   For  a  straight  TRUSS 
element  of  length  L  with  Its  center  node  at  L/2,  IP  1  Is 
located  at  0.29L,  IP  2  Is  at  L/2,  and  IP  3  Is  at  0.7IL. 
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Table  E. 1 
Normally  consolidated  soil  properties 


BBBBBCBSSSCrBEB 
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0.32 

Cc 
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Cr 
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0.04 
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0.7 

v' 

m 

0.3 
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0.47 
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s 
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Table  E.2 
Cap  parameters  for  normally  consolidated 
and  overconsol (dated  foundation  soil 

Property  Parameter   Normally        Ovei — 
number*    name    consolidated  consolidated 

BBCEBBCBBBeCEEBEEBBEBnECKE-BBEBBe-BBEBEEBEBBEBBE 


1 

Kl 

98.  0 

98.0 

2 

K2 

1.0 

1.0 

3 

AP 

2. 116  ksf 

2.1 16  ksf 

4 

Kml 

n 

1.0  ksf 

1.0  ksf 

5 

Gl 

0.0 

0.0 

6 

G2 

0.46 

0.46 

7  - 

8 

9 

10 

11 

12 

AM 

AC 

AW 

AD 

R 

XL 

(a) 
(k) 

(W) 
(D) 

(V 

0.163 

0.0001  ksf 
0.  146 
-0.996/x^ 
0.747    ° 
♦  1.0 

0.172 

0.0001  ksf 

0.146 

1.19 

0.747 

-0.74 

13 
14 
15 

TENCUT 
STATE 
Aj  <Y> 

0.1  ksf 
1.0 
0.053  kef 

0.1  ksf 
2.0 
0.053  kef 

16 

A2 

0.0 

0.0 

17 

K 
o 

0.47 

0.6 

16 
19 
20 
21 
22 

SW 1 TCH 
FAC  (6) 
TIMEON 
WGT 
BOUY 

0.0 

10.0 

0.0 

0.0 

0.0 

0.0 

10.0 

0.0 

0.0 

0.0 

•  Corresponds  to  property  number  used  In 
Appendix  D 
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Table  E.3 
Parameters  for  embankment  fill 


Property 

Parameter 

Value 

number* 

name 

SEBEEESESESE 

BBEESBBEBSSSE 

ESBBEBBBBBBBBBBBBB 

1 

K| 

190.0 

2 

K2 

0.65 

3 

A 
P 

2.1 16  ksf 

4 

mln 

20.0  ksf 

5 

G. 

0.0 

6 

G2 

0.33 

7 

AM  (o) 

0.17 

8 

AC  (O 

0.0001  ksf 

-   9 

AW  (W) 

1.0 

10 

AD  (D) 

1.0 

1  1 

R 

,'°     10 
-1.0xl0,u  ksf 

12 

XL  <t  ) 
o 

13 

TENCUT 

0.1  ksf 

14 

STATE 

2.0 

15 

Aj  (y) 

0.0 

16 

A2 

0.0 

17 

K 
o 

0.0 

18 

SWITCH 

0.0 

19 

FAC  (B) 

0.0 

20 

T1MEON 

•  • 

21 

WGT 

0.125  kef 

22 

BOUY 

0.0624  kef 

BEEEEEEEEEEE 

SEEESESSESESS 

•  Corresponds  to  property  number  used  In 

Appendix  D 
••See  section  on  selection  of  FE  mesh 


-    Ell    - 


II 


o 


o 

CO 

s 

• 

o 

X 

t- 

>- 

4-> 

■J 

_J 

X 

< 

t 

s 

p 

K 

o 

0> 

z 

+> 

eg 


I 

X 

UJ 


Id 


-  E12  - 


Program  Input: 

Input-  •m.acn,sJ2fn,rkc,abb   -then  hit  return 
Remember  that  acn  and  sJ2fn  are  negative 
0.163  0.0  -0.339  0.47  5.26 


Program  output; 

INPUT  TEST  PARAMETERS: 


note   •••  slope  •  0  error  when  computing  mex(xf/xo) 


am 

0. 16300 

acn 

0.00000 

sJ2fn 

-0.33900 

rkc 

0.47000 

abb 

5.26000 

•••  si 

ope 

B 

0  error  w 

LATED 

CAP 

PARAMETERS: 

dxo 

B 

-0.99562 

wab 

B 

2.71830 

r 

m 

0.74691 

xfxo 

m 

1.00882 

rmax 

B 

0 

.  15857E+06 

xfxomx 

B 

1. 10786 

Figure  E.2.   Example  Input  and  output  for  program  CPCALC. 


EXAMPLE  EMBANKMENT-CAP  PARAMETERS  NORMALLY  CONSOLIDATED  CLAY 

0  4 

98.0  1.0  2.116  1.0  0.0  0.46  .183  0.0001 

0.146    0.907    0.747   +1.0    0.1     1.0    0.0   -0.526   0.47 

0.0  10.0  0.0  0.0  0.0 

stress 

0.005  20  0 

0  60  0.  -0.0005  0.  0. 

0  39  0.  -0.0002  0.  0. 

0  0  0.  0.  0.  0. 


Figure  E.3.   Example  Input  for  program  CAP. 
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Figure  E.4.   J^^vo  end  Au/avo  VS*  **'*'  *tra!n 

calculated  by  program  CAP  for  normally 
consolidated  soli. 
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Figure  E.5.   Finite  element  mesh  for  example  embankment, 
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Figure  E.6  Stress-strain  curve  for  tensile  reinforcement. 


APPENDIX  E. 1 
OUTPUT  FROM  PROGRAM  CAP 
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0.3192E+00 
0.3197E+00 
0.3203E400 
0.3208E+00 
0.32J3E400 
0.3219E+00 
0.3224E+00 
0.3229E400 
0.3235E+00 
0.3240E+00 
0.3249E+00 
0.3252E+00 
0.3257E400 
0.3265E+00 
0.3268E+00 
0.3276E400 
0.3279E+00 
0.3286E+00 
0.3289E+00 
0.3297E400 
0.3303E+00 
0.3308E+00 
0.3313E+00 
0.3319E+00 
0.3324E+00 


E20    - 


OUTPUT    FROM    CAP    (CQNTI NIIED1 


CONTENTS  OF  PLOT. I  (CONTINUED): 


0. 1269E-02 
0.130IE-02 
0.133SE-02 
0.1370E-02 
0. I409E-02 
0. I442E-02 
0. 14B1E-02 
0. 1520E-02 
0. I560E-02 
0. 1606E-02 
0. 1660E-02 
0. 1709E-02 
0.1759E-02 
0.18ME-02 
0. 18B0E-02 
0. 1967E-02 
0.2093E-02 
0.2323E-02 
0.2991E-02 
0.5977E-02 
0. 1 195E-01 
0.3581E-01 
0.29B6E+00 


0.2B83E+00 
0.2887E+00 
0.2892E400 
0.2B96E+00 
0.2903E+00 
0.2906E+00 
0.2912E+00 
0.2917E+00 
0.2921E+00 
0.2925E+00 
0.2929E+00 
0.2932E+00 
0.2934E+00 
0.2936E+00 
0.2937E400 
0.2939E+00 
0.2940E+00 
0.2940E+00 
0.2941E400 
0.2941E+00 
0.2941E+00 
0.2941E+00 
0.2941E+00 


0.3224E-01 
0.3309E-01 
0.3397E-01 
0.3487E-01 
0.3579E-01 
0.3691E-01 
0.3784E-01 
0.3899E-01 
0.4015E-01 
0.4135E-01 
0.4261E-01 
0.4352E-01 
0.4434E-01 
0.4513E-01 
0.4592E-01 
0.4673E-01 
0.4758E-01 
0.4847E-01 
0.4934E-01 
0.5000E-01 
0.5042E-01 
0.5083E-01 
0.5122E-01 


0.23I6E+01 
0.2320E401 
0.2325E+01 
0.2330E+01 
0.2337E401 
0.2341E+01 
0.2347E401 
0.2353E+01 
0.2359E+01 
0.2364E+0I 
0.2369E+01 
0.2374E+01 
0.2377E40J 
0.2381E+01 
0.2384E+01 
0.2387E+01 
0.2390E+01 
0.2392E+01 
0.2393E+01 
0.2393E+01 
0.2393E+01 
0.2393E+0I 
0.2393E401 


0.3329E+00 
0.3334E400 
0.3339E400 
0.3344E+00 
0.3352E+00 
0.33S6E400 
0.3362E400 
0.3368E+00 
0.3373E+00 
0.3378E+00 
0.33B2E+00 
0.33B5E+00 
0.33B8E+00 
0.3390E+00 
0.3392E+00 
0.3393E+00 
0.3394E+00 
0.3395E+00 
0.3396E+00 
0.3396E+00 
0.3396E+00 
0.3396E+00 
0.3396E+00 


CONTENTS  OF  PLOT. 2: 


O.OOOOE+00 
0. 1102E-04 
0.3355E-04 
0.5992E-04 
0.8068E-04 
0. 1061E-03 
0. I269E-03 
0. 1527E-03 
0. I740E-03 
0.2003E-03 
0.2222E-03 
0.2491E-03 
0.2717E-03 
0.2993E-03 
0.3225E-03 
0.3445E-03 
0.3722E-03 
0.3962E-03 
0.4190E-03 
0.4417E-03 
0.4645E-03 
0.4875E-03 
0.5I09E-03 


0.3866E+00 
0.3866E+00 
0.3866E+00 
0.3865E+00 
0.3865E+00 
0.3864E+00 
0.3864E+00 
0.3863E+00 
0.3863E+00 
0.3663E+00 
0.3862E+00 
0.3862E+00 
0.3861E+00 
0.3861E+00 
0.3660E+00 
0.3860E+00 
0.38S9E+00 
0.3858E+00 
0.3857E+00 
0.3857E+00 
0.3856E+00 
0.3856E+00 
0.3855E+00 


0.1394E+00 
0.1394E+00 
0. 1396E+00 
0. 1400E+00 
0. 1401E+00 
0.1405E+00 
0. I406E+00 
0. 1410E+00 
0. 14] 1E+00 
0. 1414E+00 
0. 14J6E+00 
0. 1419E+00 
0.1421E+00 
0. 1424E+00 
0. I427E400 
0. 1428E+00 
0. 1432E+00 
0. 1434E+00 
0. 1437E+00 
0. I439E+00 
0. 1442E+00 
0. I444E+00 
0. 1447E+00 


0.1020E401 
0. 1020E+01 
0. 1020E+01 
0. 1020E+01 
0. 1019E+01 
0.1019E+01 
0. I019E+01 
0. 1018E+01 
0. 1018E+01 
0. I017E+0) 
0. 1017E+01 
0. 1017E+01 
0. 1016E+01 
0. 1016E+01 
0. 1015E+01 
0. 1015E+01 
0. 1014E+01 
0. 1014E+01 
0. 1014E+01 
0. 1013E+01 
0. I013E40I 
0. 1012E+01 
0. I012E401 


0. 1610E+00 
0. 1610E+00 
0.I6I2E+00 
0. 16I6E+00 
0. 1618E+00 
0.1622E+00 
0. 1624E400 
0.1628E+00 
0. 1629E+00 
0. 1633E+00 
0. I635E+00 
0.1639E+00 
0. 1641E+00 
0. 1645E+00 
0. 1648E+00 
0.1649E+00 
0.1653E+00 
O.I656E+00 
0.1659E+00 
0. 1662E+00 
0.I665E+00 
0. 1668E+00 
0.I671E+00 


-  E21  - 


OUTPUT  FROH  CAP  (CONTINUED) 


CONTENTS  OF  PLOT. 2  (CONTINUED): 


0.5347E-03 
0.5588E-03 
0.5634E-03 
0.6085E-03 
0.6339E-03 
0.6599E-03 
0.6864E-03 
0.7134E-03 
0.7409E-03 
0.7690E-03 
0.7977E-03 
0.8271E-03 
0.8571E-03 
0.8793E-03 
0.9062E-03 
0.9360E-03 
0.9595E-03 
0.9882E-03 
0.  1012E-02 
0.  1042E-02 
0. 1068E-02 
0. 1099E-02 
0.  1127E-02 
0.  1  154E-02 
0.  1  181E-02 
0.  1209E-02 
0.  1238E-02 
0.1269E-02 
0.I301E-02 
0.  1335E-02 
0.  1370E-02 
0.  1409E-02 
0.  1442E-02 
0.  1481E-02 
0. I520E-02 
0. 1560E-02 
0. I606E-02 
0. 1660E-02 
0.  1709E-02 
0.I759E-02 
0. I814E-02 
0. 1880E-02 
0. I967E-02 
0.2093E-02 
0.2323E-02 
0.2991E-02 
0.5977E-02 
0. 1 195E-01 
0.3581E-01 
0.2986E+00 


0.3854E+00 
0.3853E400 
0.3853E+00 
0.3852E+0D 
0.3B51E+00 
0.3850E+00 
0.3849E400 
0.3B49E+00 
0.3B48E400 
0.3847E+00 
0.3646E+00 
0.3845E+00 
0.3843E+00 
0.3842E+00 
0.3841E+00 
0.3B40E+00 
0.3839E+00 
0.3837E+00 
~0.3836E+00 
0.3834E+00 
0.3834E+00 
0.3832E+00 
0.3830E+00 
0.3828E+00 
0.3827E+00 
0.3825E+00 
0.3823E+00 
0.3821E+00 
0.3819E4-00 
0.3817E+00 
0.3815E+00 
0.381 1E  +  00 
0.3809E400 
0.3805E+00 
0.3B02E+00 
0.3798E+00 
0.3795E+00 
0.3790E+00 
0.3787E400 
0.3784E+00 
0.3781E+00 
0.3778E+00 
0.3775E400 
0.3772E+00 
0.3769E+00 
0.3768E+00 
0.3768E+00 
0.3768E+00 
0.3767E+00 
0.3768E+00 


0.  I449E400 
0.1452E+00 
0.  14S4E+00 
0. 1456E+00 
0. 1459E*00 
0.  1461E  +  00 
0. I464E+00 
0. 1466E+00 
0. 1469E+00 
0.  1471E+00 
0. 1473E400 
0. I476E400 
0. 1480E+00 
0.  1481E  +  00 
0. 1483E+00 
0. 1487E+00 
0. 1489E+00 
0. 1492E+00 
0. 1494E+00 
0. 1497E+00 
0. 149BE+00 
0. 1502E400 
0. 1504E+00 
0.1507E+00 
0.  I509E400 
0.1512E+00 
0. 1514E+00 
0.1516E+00 
0. 1519E+00 
0.1521E+00 
0. 1523E+00 
0.1527E+00 
0. 1529E+00 
0.1532E+00 
0.  1534E+00 
0.1536E+00 
0. 1539E+00 
0.1541E+00 
0. 1542E400 
0.  1543E+00 
0. 1544E+00 
0. 1545E+00 
0. I546E+00 
0. 1546E+00 
0. 1547E400 
0.  1547E+00 
0.  1547E+00 
0. I547E+00 
0. 1547E+00 
0. I547E+00 


0. 101 1E+01 

0. 101 1E401 

0.  10I0E  +  01 

0.1010E+01 

0.  1009E401 

0. 1009E+01 

0. 1008E+01 

0. 1008E+01 

0. 1007E+01 

0.  1007E+01 

0. 1006E+01 

0.  1006E+01 

0. 1005E+01 

0.1005E+01 

0. 1004E+01 

0. 1003E+01 

0. 1003E+01 

0.1002E+01 

0. 1002E+01 

0. 1001E+01 

0.  1000E+01 

0.9993E+00 

0.9985E+00 

0.9978E+00 

0.9970E400 

0.9963E+00 

0.9955E+00 

0.9947E+00 

0.9939E+00 

0.9930E+00 

0.9921E+00 

0.9906E+00 

0.9899E+00 

0.9884E+00 

0.9871E+00 

0.9859E+00 

0.9845E+00 

0.9831E+00 

0.9819E+00 

0.9809E+00 

0.9799E+00 

0.9789E+00 

0.9779E+00 

0.9769E+00 

0.9760E+00 

0.9756E+00 

0.9756E400 

0.9756E+00 

0.9756E400 

0.9756E+00 


O.I673E*00 
0.I676E+00 
0, 1679E+00 
O.I682E+00 
0.1685E+00 
0.1687E*00 
0.1690E+00 
0.1693E+00 
0. 1696E+00 
0. 1699E+00 
0.1701E+00 
0.1704E+00 
0.1709E+00 
0.171 1E+00 
0.1713E+00 
0.1717E+00 
0. 1719E+00 
0.1723E+00 
0.1725E+00 
0.1729E+00 
0.1730E+00 
0.1734E+00 
0.1737E+00 
0.1740E+00 
0.1743E+00 
0.1746E+00 
0.1748E+00 
0.1751E+00 
0.1754E+00 
0.1756E+00 
0.1759E+00 
0.1763E+00 
0.1765E+00 
0.1769E+00 
0.1771E400 
0.1774E+00 
0.1777E+00 
0.1779E+00 
0.1781E400 
0.1782E+00 
0.1783E+00 
0.  1764E  +  00 
0.1785E^00 
O.I785E+00 
0.1786E+00 
0.1786E+00 
0.1786E*00 
0.1786E+00 
0.1786E+00 
0. 1786E+00 


APPENDIX  E.2 
INPUT  TO  PROGRAM  AUTOGEN 
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INPUT  TO  AUTOGEN 

Items  In  parentheses  are  comments  added  to  explain  the 
commands  and  are  not  part  of  the  Input  file. 


193 

52 

(card  1,  number  of  nodes  and  elements) 

NFAP 

(control  Information) 

EXAMPLE 

EMB 

-  90*  BASE  WIDTH  -  2:1  SLOPE  -  WEAK  CRUST  - 

1.875'  LIFTS 

193 

0 

2    I 

1    2 

1     15930    1 

1 

0 

1 

15    1 

.01 

0. 

•  i 

5 

1. 

0 

1 

0 

JOINT 

(Input  coordinates  of 

selected  nodes) 

1 

0. 

0. 

10 

7.5 

0. 

136 

0. 

0.9375 

148 

7.5 

1.875 

184 

35.625 

1.875 

187 

43.125 

0.9375 

MSHGN 

(generate  foundation  and  some  embankment 

nodes) 

1 

10 

14 

15. 

9    1 

-3.75 

10 

9 

14 

15. 

5    1 

-7.5 

136 

3 

18 

15. 

12    1 

0.9375 

148 

2 

18 

15. 

6    1 

1.875 

JTGEN 

(generate  more  embankment  nodes) 

179 

5 

1 

-1.675 

0.9375 

184 

3 

1 

-1.875 

1.875 

187 

7 

1 

-1.875 

0.9375 

JOINT 

(redefine  nodes  170  end  171) 

170 

20.625 

9.375 

171 

18.75 

11.25 

JTFIX 

(specify  fl 

xed  degrees 

»  of  freedom) 

0 

1 

0 

1    6 

127   134 

136   147 

1 

1 

1 

9   14 

23   28 

37   42   51   56 

65   70   79   64 

1 

1 

1 

98  107 

112   121 

126  135 

NFAP 

(load  control  and  2-d 

element  data) 

135 

7 

4 

0 

2  52 
■ 

4 

8 

2 

10  10  22 

12 

98. 

1. 

.0 

2.116 

1. 

0.0       0.46 

0.172 

0.0001 

0.146 

1. 

.19 

0.747 

-0.740 

0.1       2.0 

0.053 

0.0 

0.6 
2 
96. 

0. 

t 

-10. 

0. 

0.       0. 

1. 

.0 

2.116 

1. 

0.0       0.46 

0.163 

0.0001 

0.146 

0. 

.606 

0.747 

4J. 

0.1       1. 

0.053 

0.0 

0.47 
3 

96.0 

0. 

10. 

0. 

0.        0. 

1. 

.0 

2.116 

1. 

0.0      0.46 

0.183 

0.0001 

0.146 

0. 

.484 

0.747 

♦1.0 

0.1       1.0 

0.053 

0.0 

0.47 

4 
98.0 

0. 

.0 

10. 

0.0 

0.0      0.0 

0.0 

1. 

.0 

2.116 

1. 

0.0      0.46 

0.183 

0.0001 

0.146 

0 

.346 

0.747 

♦  1.0 

0.1      1.0 

0.053 

0.0 

0.47 

0. 

.0 

10. 

0.0 

0.0      0.0 

0.0 

93 
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INPUT    TO   AUTQGEN    (CONTINUED) 


5 

190. 

0 

.65 

2 

116 

20 

• 

0.0 

0.33 

1. 

i 

ft 

1 

t 

0. 

m 

01 

1.E9  .01 
.125 

2. 

.0624 

6 

190. 

0 

.65 

2. 

116 

20 

. 

0.0 

0.33 

1. 

1 

1 

1. 

1.E9  .01 
.125 

2. 
.0624 

7 

190. 

0 

.65 

2. 

116 

20 

• 

0.0 

0.33 

1. 

1 

• 

1 

i 

2. 

m 

I.E9  .01 
.125 

2. 

.0624 

8 

190. 

0 

.65 

2 

116 

20 

, 

0.0 

0.33 

1. 

1 

ft 

1 

3 

• 

1.E9  .01 
.125 

2. 

.0624 

9 

190. 

0 

.65 

2. 

116 

20 

, 

0.0 

0.33 

1. 

1 

• 

1 

4 

I.E9  .01 
.125 

2. 

.0624 

10 

190. 

0 

.65 

2. 

116 

20 

a 

0.0 

0.33 

1. 

1 

» 

1 

5 

1.E9  .01 
.125 

2. 

.0624 

ELEME 

(define  selected 

2-d 

elements) 

1 

1 

15 

3 

17 

10 

2 

11 

16 

2 

2 

17 

5 

19 

11 

4 

12 

18 

3 

3 

19 

7 

21 

12 

6 

13 

20 

4 

4 

21 

9 

23 

13 

8 

14 

22 

37 

5 

155 

i: 

1 

15 

148 

136 

10 

154 

38 

6 

157 

1: 

137 

155 

149 

138 

148 

156 

39 

7 

159 

h 

139 

157 

150 

140 

149 

158 

40 

8 

161 

u 

141 

159 

151 

142 

150 

160 

41 

9 

163 

u 

143 

161 

152 

144 

151 

162 

42 

10 

165 

u 

145 

163 

153 

146 

152 

164 

43 

5 

173 

1! 

15 

29 

166 

154 

24 

172 

49 

5 

186 

r 

29 

43 

184 

172 

38 

187 

50 

6 

190 

p 

173 

188 

185 

174 

184 

189 

51 

7 

192 

r 

175 

190 

186 

176 

185 

191 

0.17 


0.17 


0.17 


0.17 


0.17 


0.17 


0.0001 


0.0001 


0.000! 


0.0001 


0.0001 


0.0001 
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INPUT  TO  AUTQGEN  (CONTINUED I 


52 
179 

ELSTR 

I 
2 
3 

4 
ELDUP 

38 
NFAP 
1 
I 


-1 


0. 
0. 
0. 
0. 
0. 
0. 

0. 

LOAD 

FOLL 


I 

2 
3 

1 

2 
3 
4 
5 
6 
7 


i   e 

179  177   192   179  17B  1B6   193 

(generate  foundation  elements) 
9    4    14    14    14    14    14    14    14    14 
9    4    14    14    14    14    14    |4    |4    14 
9    4    14    14    14    14    14    14    14    14 
9    4    14    14    14    14    14    14    14    14 

(generate  embankment  elements  44  through  48) 
42    6   18 

(truss  element  data  and  load  curves) 


15 

29 

2 


-0.0001 
10 
24 

38 

-I. 

0. 
0. 
0. 
0. 
0. 
0. 


0. 

0. 


15 
29 

43 


10. 

I. 

1. 

2. 

3. 

4. 


-1. 
-I. 

0. 
0. 
0. 
0. 


0. 

.0001 


(load  curve 
( load  curve 

10. 

2. 

3. 

4. 

5. 


60. 


1) 
I) 


10. 
10. 
10. 
10. 
10. 


5.        0.       6. 

(generate  nodal  loads  to  balance  Initial  stress  In  foundation) 
1        1 
1       36    0.053         2        3 

(last  card) 


APPENDIX  E.3 
INPUT  TO  PROGRAM  PS-NFAP 
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INPUT    TO    PS-NFAP 


tXAHPLE 

EHB  • 

•  90' 

BASE  WIDTH  -  2:1  SLOPE  -  WEAK 

CRUST  - 

1.875'  LI 

193 

0 

2 

1    1    2 

1     1S930    1 

1 

0 

1 

IS 

1 

.01 

0. 

.5 

1 

• 

0 

1 

0 

ID 

0 

0 

(note:  D  added 

0.000 

0.000 

b.ooo 

2 

0 

0 

with  text  editor 

;  0.000 

0.000 

-3.750 

3 

0 

0 

not  generated 

0.000 

0.000 

-7.500 

4 

0 

0 

by  AUTOGEN) 

0.000 

0.000 

-11.250 

5 

0 

0 

0.000 

0.000 

-15.000 

6 

0 

0 

0.000 

0.000 

-18.750 

7 

0 

0 

0.000 

0.000 

-22.500 

e 

0 

0 

0.000 

0.000 

-26.250 

9 

0 

1 

0.000 

0.000 

-30.000 

10 

0 

0 

0 

0.000 

7.500 

0.000 

11 

0 

0 

0 

0.000 

7.500 

-7.500 

12 

0 

0 

0 

0.000 

7.500 

-15.000 

13 

0 

0 

0 

0.000 

7.500 

-22.500 

14 

0 

1 

1 

0.000 

7.500 
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JNPUT 

TO 

PS-NFAP     (CONTINUED) 

1 

3 

I-.4969E+00 

3 

I0.I987E*0I 

3 

1-.9937E+00 

3 

I0.I967E+01 

3 

I-.9937E+00 

3 

10.1987E401 

3 

1-.993BU00 

3 

10. 1967E+01 

3 

I-.4969E+00 

10 

3 

10.1987E+01 

11 

3 

I0.3975E+01 

12 

3 

I0.3975E+0I 

13 

3 

I0.3975E+01 

14 

3 

10. 1987E+01 

15 

3 

I-.9937E400 

16 

3 

10.3975E+0I 

17 

3 

1-.1987E+01 

18 

3 

I0.3975E+01 

19 

3 

1-.1987E+01 

20 

3 

10.3975E+0I 

21 

3 

1-.1987E+01 

22 

3 

J0.3975E+01 

23 

3 

1-.9937E+00 

24 

3 

10. 1987E+01 

25 

3 

10.3975E+01 

26 

3 

10.3975E+0I 

27 

3 

10.3975E+01 

28 

3 

10. 1987E+01 

29 

3 

1-.9937E+00 

30 

3 

I0.3975E+01 

31 

3 

1-.1987E+01 

32 

3 

I0.3975E+01 

33 

3 

I-.1987E+01 

34 

3 

I0.3975E+0I 

35 

3 

1-.1987E+01 

36 

3 

I0.3975E+0J 

37 

3 

1-.9937E+00 

38 

3 

10. 1987E+01 

39 

3 

10.3975E+01 

40 

3 

I0.3975E+0I 

41 

3 

I0.3975E+01 

42 

3 

I0.I987E+0I 

43 

3 

1-.9938E+00 

44 

3 

I0.3975E40I 

45 

3 

1-.I987E+01 

46 

3 

10.3975E+01 

47 

3 

I-.I987E+01 

48 

3 

I0.3975E*OI 

49 

3 

1-.I987E+0I 

50 

3 

10.3975E401 

51 

3 

1-.9937E+00 

52 

3 

10.1987E+0I 

53 

3 

I0.3975E+OI 

-    E33   - 


INPUT 

TO  PS-NFAP  (CONTINUED) 

54 

3 

10.3975E+01 

55 

3 

I0.3975E*01 

56 

3 

10. 1987E+01 

57 

3 

I-.9937E+00 

58 

3 

10.3975E+01 

59 

3 

I-.I987E+01 

60 

3 

10.3975E+01 

61 

3 

1-.1987E+01 

62 

3 

I0.3975E+01 

63 

3 

I-.I987E+0I 

64 

3 

I0.3975E+01 

65 

3 

1-.9937E+00 

66 

3 

10.1987E+01 

67 

3 

10.3975E+01 

68 

3 

I0.3975E+01 

69 

3 

10.3975E401 

70 

3 

10.1987E+01 

71 

3 

I-.9938E+00 

72 

3 

I0.3975E+01 

73 

3 

I-.1987E+01 

74 

3 

10.3975E+01 

75 

3 

I-.19B7E+01 

76 

3 

10.3975E+01 

77 

3 

1-.1987E+0I 

78 

3 

10.3975E+01 

79 

3 

1-.9937E+00 

80 

3 

10.19B7E+01 

81 

3 

I0.3975E401 

82 

3 

10.3975E+01 

83 

3 

I0.3975E40I 

84 

3 

10.I987E+01 

85 

3 

1-.9937E400 

86 

3 

10.3975E+01 

87 

3 

1-.1988E401 

88 

3 

I0.3975E+0! 

89 

3 

1-.1987E+01 

90 

3 

10.3975E+01 

91 

3 

I-.I987E+0I 

92 

3 

10.3975E401 

93 

3 

I-.9937E400 

94 

3 

10. 1988E+01 

95 

3 

10.3975E401 

96 

3 

I0.3975E401 

§7 

3 

10.397SE401 

96 

3 

10.1988E401 

99 

3 

I-.9937E400 

too 

3 

10.3975E+01 

101 

3 

1-.1988E+0I 

102 

3 

I0.3975E40I 

103 

3 

I-.1987E40I 

104 

3 

I0.3975E401 

105 

3 

I-.I988E401 

106 

3 

I0.3975E401 

E34 


INPUT 

TO   PS 

-NFAP    (CONTINUED) 

107 

3 

1-.9938E+00 

108 

3 

10. 1967E+0I 

109 

3 

I0.3975E+0I 

110 

3 

I0.3975E*01 

111 

3 

I0.3975E+01 

112 

3 

10.I987E+0) 

113 

3 

)-.993BE*00 

114 

3 

10.3975E401 

115 

3 

I-.I987E+0I 

116 

3 

I0.3975E401 

117 

3 

1-.1987E+01 

118 

3 

I0.3975E+01 

119 

3 

I-.1987E+01 

120 

3 

I0.3975E401 

121 

3 

1-.9937E+00 

122 

3 

I0.1987E+01 

123 

3 

10.3975E40! 

124 

3 

I0.3975E+01 

125 

3 

10.3975E+01 

126 

3 

10.I988E+01 

127 

3 

1-.4969E+00 

128 

3 

10.1987E+01 

129 

3 

1-.9937E+00 

130 

3 

I0.1987E+0I 

131 

3 

1-.9937E+00 

132 

3 

I0.I987E+01 

133 

3 

1-.9937E+00 

134 

3 

10.1987E+0I 

135 

3 

1-.4969E+00            TOT/ 

Z  -  DIR  LOAD  -   0.2146500E+03  NCARD  -      135 


APPENDIX  E.4 
OUTPUT  FROM  PROGRAM  NFMINX 
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APPENDIX  E.5 
OUTPUT  FROM  PROGRAM  PS-NFAP 
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